DY T7010 « 


The Weizmann Science Press of Israel e Jerusalem p°>yi7> © Sxawoa MIIdIDDA VARA YTHI OMoID> Y 


Volume 10C, Number 1-2, June 1961 


BULLETIN 


OF THE RESEARCH COUNCIL 
OF ISRAEL 


Section ¢€ 
TECHNGLOGY 


Bull. Res. Counc. of Israel. C. Techn. 


Incorporating the Scientific Publications of the 
Technion — Israel Institute of Technology, Haifa 


A century of flow through porous media and its unsolved problem 


S. Irmay 
Cross stresses in laminar flow of gases in accordance with Maxwell's 
dynamical theory M. Reiner 
A study of stability at supercritical speeds H. Wolski 


Energy levels in modes of free extensional vibrations of a finite plate 

A. Foux and N. Davids 

On the dynamical derivative of physical tensors Z. Karni 

Utilisation of digital computers in structural analysis J. Loewinger 
Stresses and deformation in internally loaded containers 

H. Ilberg and Y. Yarnitsky 

Beams of varying cross section in bending M. Mitelman 

Elastic analysis of composite walls—A general theory S. Rosenhaupt 

A numerical method for the determination of stresses in walls with 


rectangular openings G. Miiller 
Strength and durability of mortars made with an admixture of local 
pozzuolana R. Malinowski 


The recovery of crimp wool fibres /. Banbaji, E. Alexander and M.Lewin 
A new method for measuring and defining crimp of textile fibres 

M. Shiloh, D. Mejzler and E. Alexander 
The stability of the motion of a fluid confined between two concentric 


cylinders rotating in opposite directions L. Rintel 
Characteristics of turbulent flow with suspensions of solid particles 

C. Elata 
Some experimental investigations on natural frequencies of ring 
formed bodies M. Ehrenreich and S. Braun 
The vibrations of a membrane in contact with air M. Bentwich 
On curved tubes with variable cross section A. Kornecki 
Fracture of solids under combined stresses E. Ben-Zyvi 


Shear strength analysis of specimens having a high diameter-height 
ratio E. Shklarsky and M. Livreh 


Digitized by the Internet Archive 
in 2023 


https://archive.org/details/israel-journal-of-technology 1961-06 1 Oc_1 -2 


a re 


The Weizmann Science Press of Israe! e Jerusalem o°%w17> © Sxawon ADVI VIVA VINA ODIO? 131 TON 


Volume 10C, Number 1-2, June 1961 


BULLETIN 


OF THE RESEARCH COUNCIL 
OF ISRAEL 


Section ¢€ 
TECHNOLOGY 


Bull. Res. Counc. of Israel. C. Techn. 


Incorporating the Scientific Publications of the 
Technion — Israel Institute of Technology, Haifa 


A century of flow through porous media and its unsolved problem 


S. Irmay 
Cross stresses in laminar flow of gases in accordance with Maxwell’s 
dynamical theory M. Reiner 
A study of stability at supercritical speeds H. Wolski 


Energy levels in modes of free extensional vibrations of a finite plate 

A. Foux and N. Davids 

On the dynamical derivative of physical tensors Z. Karni 

Utilisation of digital computers in structural analysis J. Loewinger 
Stresses and deformation in internally loaded containers 

H. Ilberg and Y. Yarnitskn 

Beams of varying cross section in bending M. Mitelmay 

Elastic analysis of composite walls—A general theory S. Rosenhaupt 

A numerical method for the determination of stresses in walls with 


rectangular openings G. Miiller 
Strength and durability of mortars made with an admixture of local 
pozzuolana R. Malinowski 


The recovery of crimp wool fibres. Banbaji, E. Alexander and M.Lewin 
A new method for measuring and defining crimp of textile fibres 

M. Shiloh, D. Mejzler and E. Alexander 
The stability of the motion of a fluid confined between two concentric 


cylinders rotating in opposite directions L. Rintel 
Characteristics of turbulent flow with suspensions of solid particles 

C. Elata 
Some experimental investigations on natural frequencies of ring 
formed bodies M. Ehrenreich and S. Braun 
The vibrations of a membrane in contact with air M. Bentwich 
On curved tubes with variable cross section A. Kornecki 
Fracture of solids under combined stresses E. Ben-Zvi 


Shear strength analysis of specimens having a high diameter-height 
ratio E. Shklarsky and M. Livneh 


BULLETIN 
OF THE RESEARCH COUNCIL 
OF ISRAEL 


MIRIAM BALABAN 
Editor 


EDITORIAL BOARDS 


SECTION A 
CHEMISTRY 


Y. AVIDOR 
E. D. BERGMANN 
H. BERNSTEIN 
M,. R. BLOCH 
E, KATCHALSKI 
A. KATZIR (KATCHALSKY) 
G. STEIN 
(Chairman, 
Tsrael Chemical Society) 


SECTION B 
ZOOLOGY 


H. MENDELSSOHN 
K. REICH 
L. SACHS 
A. YASHOUV 


SECTION C 
TECHNCLOGY 


A. BANIEL 
J. BRAVERMAN 
A. DE LEEUW 
M, LEWIN 
M. REINER 
A. TALMI 
E. GOLDBERG, Technion 
Publications Language Editor 


SECTION D 
BOTANY 


N. FEINBRUN 
N. LANDAU 
H. OPPENHEIMER 

T. RAYSS 
J, REICHERT 
M. ZOHARY 


SECTION E 
EXPERIMENTAL MEDICINE 


S. ADLER 
A. DE VRIES 
A, FEIGENBAUM 
M. RACHMILEWITZ 
B. ZONDEK 


SECTION F 


MATHEMATICS and PHYSICS 


A. DVORETZKY 
J. GILLIS 

F, OLLENDORFP 
G. RACAH 


SECTION G 
GEO-SCIENCES 


G. DESSAU 
J. NEUMANN 
L. PICARD 


NOTICE TO. CONTRIBUTORS 


Contributors to the Bulletin of the Research Council of Israel should conform 
to the following recommendations of the editors of this journal in preparing 
manuscripts for the press. 

Contributions must be original and should not have been published 
previously. When a paper has been accepted for publication, the author(s) 
may not publish it elsewhere unless permission is received from the Editor 
of this journal. 

Papers may be submitted in English and in French. 


MANUSCRIPT 
General 
Papers should be written as concisely as possible. MSS should be typewritten 
on one side only and double-spaced, with side margins not less than 2.5 cm 
wide. Pages, including those containing illustrations, references or tables, 
should be numbered. 
The Editor reserves the right to return a MS to the author for retyping or 
any alterations. Authors should retain copies of their MS. 


Spelling 
Spelling should be based on the Oxford Dictionary and should be consistent 
throughout the paper. Geographic and proper names in particular should be 
checked for approved forms of spelling or transliteration. 


Indications 


Greek letters should be indicated in a legend preceding the MS, as well as by a 
pencil note in the margin on first appearance in the text. 

When there is any room for confusion of symbols, they should be carefully 
differentiated, e.g. the letter ‘‘1”’ and the figure “‘1’’; ““O”’ and “0”. 


Abbreviations 
Titles of journals should be abbreviated according to the World List of Scientific 
Periodicals. 
Abstract 
Every paper must be accompanied by a brief but comprehensive abstract. 
Although the length of the abstract is left to the discretion of the author, 3 % of 
the total length of the paper is suggested. 


References 
In Sections A and C, and in Letters to the Editor in all Sections, references 
are to be cited in the text by number, e.g. ... Taylor’ ..., and are to be 


arranged in the order of appearance. 

In Sections B, D, E, and G, the references are to be cited in the text by the 
autho1’s name and date of publication in parentheses, e.g. ... (Taylor 1932)... 
If the author’s name is already mentioned in the text, then the year only 
appears in the parentheses, e.g. ... found by Taylor (1932)... The references 
in these Sections are to be arranged in alphabetical order. 

In Section F, references are to be cited in the text by number in square 
brackets, e.g. ... Taylor[3].., and are to be arranged in alphabetical order. 

The following form should be used: 

3. Taytor, G.I., 1932, Proc. roy. Soc., A138, 41. 

Book references should be prepared according to the following form: 

4. Jackson, F., 1930, Thermodynamics, 4th ed., Wiley, New York. 


TYPOGRAPHY 


In all matters of typography the form adopted in this issue should be followed. 
Particular attention should be given to position of symbols, headings, etc. 
and type specification. 


ILLUSTRATIONS 


Illustrations should be sent in a state suitable for direct photographic repro- 
duction. Line drawings should be drawn in large scale with India ink on white 
drawing paper, bristol board, tracing paper, blue linen, or blue-lined grapk 
paper. If the lettering cannot be drawn neatly by the author, he should indicate 
itin pencil for the guidance of the draftsman. Possible photographic reduction 
should be carefully considered when lettering and in other details. 

Half-tone photographs should be on glossy contrast paper. 

Illustrations should be mounted on separate sheets of paper on which the 
caption and figure number is typed. Each drawing and photograph should be 
identified on the back with the author’s name and figure number. 

sng place in which the figure is to appear should be indicated on the margin 
of the ' 


PROOFS 


Authors making revisions in proofs will be required to bear the costs thereof- 
Proofs should be returned to the Editor within 24 hours, otherwise no respon. 
sibility is assumed for the corrections of the author. 


REPRINTS 


Reprints may be ordered at the time the proof is returned. A table designating 
the cost of reprints may be obtained on request. 


Orders in America should be addressed to the Weizmann Science Press, P.O.B. 801 Jerusalem or through booksellers» 

and in England and Europe to Wm. Dawson and Sons, Ltd. Cannon House, Macklin Street, London W.C. 2, directly 

or through booksellers. Annual subscription per section (four issues): IL. 6.000 ($6.00, £2.02). Single copy: IL. 1.500 

($1.50. 12s.) — Manuscripts should be addressed: The Editor. The Weizmann Science Press of Israel, P.O.B. 801 
Jerusalem, 33 King Geotge Ave. Telephone 27844, 26345. 


THE ISRAEL SOCIETY FOR THEORETICAL AND APPLIED MECHANICS 


The Israel Society for Theoretical and Applied Mechanics was founded in 1950 
on the initiative of Professor M. Reiner. 

The purpose of the Society is to advance the development of this branch of science 
in Israel, by the concentration of scientific resources by mutual aid in research, 
and by organizing scientific gatherings on the subjects of Theoretical and Applied 
Mechanics. The Society also represents Israel in the Internation Union of Theoretical 
and Applied Mechanics (I.U.T.A.M.). 


MEMBERS OF THE EXECUTIVE COMMITTEE 


President — Prof. M. REINER 
Secretary — Myr. A. Foux 
Member —— Dr. D. ABIR 


MEMBERS OF THE PROGRAM COMMITTEE OF 
THE 10TH CONFERENCE 
Dr. D. ABIR 
Dr. A. BETSER 
Dr. J. GLUCKLICH 
Dr. Z. KARNI 
Mr. A. Foux 


PROGRAM OF THE TENTH CONFERENCE 
under the auspices of the 
TECHNION-ISRAEL INSTITUTE OF TECHNOLOGY, DIVISION OF MECHANICS 
HAIFA 
June 28-29, 1961 
28.6.1961 
09.00-09.30: Registration 


09.30-10.00: FIRST SESSION—GENERAL MEETING 
Chairman: Prof. M. REINER 
1. President’s report 
2. Secretary’s report 
3. Elections 


10.00-11.10: SECOND SESSION—MAIN LECTURE 
Chairman: Prof. M. REINER 
A Century of flow through porous media and S. IRMAY 
its unsolved problem 


11.10-13.10: 


14.15-17.35: 


29.6.1961 


THIRD SESSION—DYNAMICS OF FLUIDS 


Chairman: Dr. M. HANIN 
Cross stresses in laminar flow of gases in ac- 
cordance with Maxwell’s dynamical theory 


The stability of the motion of a fluid confined 
between two concentric cylinders rotating in 
opposite directions 


Characteristics of turbulentflow with suspensions 
of solid particles 


FOURTH SESSION—DYNAMICS OF SOLIDS 


Chairman: Prof. E. JABOTINSKY 


A study of stability at supercritical speeds 


Energy level in modes offree extensional vibrations 
of a finite plate 


Some experimental investigations on natural fre- 
quencies of ringed formed bodies 


The vibrations of a membrane in contact with air 
On the dynamical derivative of physical tensors 


M. REINER 


L. RINTEL 


H. ELATA 


H. Wo.Lsky 
A. Foux AND N. DAvips 


M. EHRENREICH AND S. BRAUN 


M. BENTWICH 
Z. KARNI 


09.30-11.30: FIFTH SESSION — PROPERTIES OF MATERIALS 


11.30-12.50: 


14.00-18.00: 


Chairman: Prof. M. REINER 
The recovery of crimp of wool fibres 


Fracture of solids under combined stresses 
Strength and durability of mortars made with an 
admixture of lccal pozzoulana 


SIXTH SESSION—COMPUTING METHODS 


Chairman: Dr. D. ABIR 
A new method for measuring and defining crimp 
of textile fibres 


Utilisation of digital computers in structural 
analysis 


J. BENBANI, E. ALEXANDER AND 
M. LEWIN 


E. BEN-ZvVI 
R. MALINOVSKI 


M. SHILoH, D. MEJZLER AND 
E. ALEXANDER 


J. LOEWINGER 


SEVENTH SESSION—ELASTICITY AND STRENGTH OF MATERIALS 


Chairman: Dr. J. GLUCKLICH 
On curved tubes with variable cross section 


Stresses and deformation in internally loaded 
containers 


Shear strength analysis of specimens having a high 
diameter-height ratio 


Beams of varying cross section in bending 


Elastic analysis of composite walls—A general 
theory 


A numerical method for the determination of 
stresses in walls with rectangular openings 


A. KORNECKI 
H.ILBERG AND Y. YARNITSKY 


E. SHKLARSKY AND M. LIvNEH 


M. MITELMAN 
S. ROSENHAUPT 


G. MUELLER 


A CENTURY OF FLOW THROUGH POROUS MEDIA AND ITS 
UNSOLVED PROBLEM 


S. IRMAY 


Hydraulic Laboratory, Department of Civil Engineering, 
Technion — Israel Institute of Technology, Haifa 


ABSTRACT 


Flow of fluids through porous media has been investigated for overacentury, yet 
many problems remain unsolved. In spite of the hundreds of papers published 
yearly on the subject, no comprehensive textbook is available, except on limited 
portions. This is due to the intricacy of such flows. 

1. The porous medium may be chemically and physically active orinert; stable 
or unstable; with constant or variable porosity; homogenous or heterogenous; 
isotropic or not. 

2. The liquid may be homogeneous or consist of several immiscible or miscible 
liquids interspersed with isolated gas bubbles or having a continuous gaseous phase 
with variable amounts of vapour. The bond of the liquid to the soil may be chem- 
ical, physico-chemical (hydratation water) or physico-mechanical (capillarity). The 
forces involved are not well known and well defined; such practical concepts as 
capillary potential need re-defining. 

3. The flow may be laminar with negligible inertial terms (Darcy flow) or important 
inertial terms (Dupuit-Forchheimer flow), or even turbulent. Little is known about 
the non-Darcy flows. 

4. Confined saturated flow can at least be studied using known mathematical 
tools ,such as theory of functions, conformal mapping and analogies. Unconfined 
flow with free boundaries constitutes a non-linear problem with an undefined bounda- 
ry. Only few solutions are known. This becomes highly complicated in the case of 
two liquids in contact (oil or sea-water and water). 

5. Unsteady flow is the field with the least reliable knowledge. Besides the 
mathematical intricacy of the complex non-linear problem with undefined boundaries 
variable in time, there is a serious doubt about the character of the flowitself, which 
seems to be unsaturated with the possibility of hysteresic phenomena on wetting and 
drying. Thermodynamic, physico-chemical and even biological factors have a highly 
disturbing effect. 

6. Some new basic solutions are presented, which may serve as a first approxi 
mation to a number of important applications. 
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CROSS-STRESSES IN THE LAMINAR FLOW OF GASES IN ACCORDANCE 
WITH MAXWELL’S DYNAMICAL* THEORY 


M. REINER 
Division of Mechanics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


In a number of papers Popper and Reiner 5:6.8:10 have described various instruments 
which show that “‘cross-stresses” are present in the laminar flow of air, i.e. stresses acting 
in the direction of the velocity of flow and, normal to it, in the direction of the velocity 
gradient. The theoretical explanation of the effect is still an open problem. The current 
kinetic theories which take into account second order effects, such as those of Chapman 
and Cowling! and of Ikenberry and Truesdell? predict results of opposite sign. i.e. pressure 
in the direction of the flow, or an increase of pressure in this direction and therefore a 
centrifugal effect and, accompanying this, a reduction of the pressure in the direction 
of the velocity gradient or, if one may say so, a negative thrust bearing. It is therefore 
interesting to note that cross-stresses such as mentioned above, resulting in centripetal 
pumping action and a self supporting thrust bearing action, seem to follow from the “‘dy- 
namical’’ theory as developed by Maxwell3. To show this, is the purpose of the present 
paper where I shall proceed by elementary methods on simplifying assumptions. While 
such procedure cannot be expected to be quantitatively correct, it should not affect 
the generality of the results in their qualitative aspects. 


We start from a theorem expressed by Maxwell? in the following words: 

‘‘Suppose the molecules to be confined in a rectangular vessel with perfectly 
elastic sides, and that they have no action on one another, so that they never strike 
one another, or cause each other to deviate from their rectilinear paths... Now 
suppose the pressure at first equal in the three directions perpendicular to the sides 
and let the dimensions a, b, c of the vessel be altered by small quantities da, 5b, dc. 

Then, if the original pressure in the direction a was p, it will become** 


5p_ _ 6a (tab) ¥ 


- This is how Maxwell himself called his theory, which we now name “‘kinetic’’. 
** It is assumed to be “at first’? equal to p also in the directions b and c, i.e. we start from a steady 
state. 


Recieved May 15, 1961. 
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Maxwell does not give any proof or explanation, apparently considering what he 
says as obvious. However, in order to understand the limitations of his theorem, 
due to the assumptions,made by him, and to estimate the amendments which have 
to be made if these assumptions do not hold, it will be useful to provide a detailed 
argument. 

In accordance with the elementary kinetic theory we consider one molecule of 
mass m moving in the direction a with velocity u. By impact with one of the sides 
normal to a when at rest, it transfers the momentum mu —(—mu) = 2mu. It returns 
to the same side after a lapse of the time 2a/u, or at a frequency u/2a, and the impact 
per unit time, or the pressure force on the side of area bc is accordingly 


2mu——s mu? 
be Daa ma 2alu 4 a (1.2) 
Therefore writing 
V = abe (1.3) 


for the volume of the vessel, we have from one molecule 


Poa = mu?/V (1.4) 
and from N molecules 
an = 
Paa = mee =nmu? (1.5) 
where 
n= af (1.6) 


is the number-density of the molecules, and # is the mean value of the velocities. 
In the steady state, the motions of the molecules have no preference for particular 
directions in space, and miu” = mv” = mw”, where v and w are the velocities in the 


directions b and c. We then have for the pressure 
2 


Paa = Pow = Pee = | La tee (1.7) 
where C, the average total velocity, is defined by 
C?im = Lm(u? + v* + w’) (1.8) 


Now, let us move one of the sides normal to a a distance da at the velocity U, which 
is smaller than u. Due to the “‘perfect elasticity’ of the side, the velocity of the 
molecule is reversed by collision. Relative to the moving side, the velocity of the 
molecule is u — U before collision, and —(u—- U) = —u+U after collision. Re- 
lative to a fixed system it is therefore 


(—u+U)+U=2U—-u (1.9) 
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and by one collision the velocity with which the molecule hits the wall is reduced 
by 2U. This process goes on for the interval T, where 


TU = 6a (1.10) 


During this interval the molecule hits the side uT/2a times = 6a/2aU. After the 
total expansion, the molecule therefore hits the wall with the velocity 


da da 
u—2U > =u(1-—) (1.11) 
The pressure is now in analogy with (1.2) 
da 
2 —— ee 
ae ae eee (=) 1.12 
Pas = Ai a) ( ) 
fen 
or to the first order 
s bas L6bMu OE 
Peo Tuh-Als Serenieae a) 
3 da 6b dc 
By = p (1-—-3 [-=) (1:13) 
S da ob dc 
Pax rep. ane yep) 
where 
p = nmi? = nmv? = nmw? 
in accordance with Eq. (1.5) 
N 
ote that Sa 6b bc _ 8V 
he tia, oars (1-15) 


is equal to the cubical dilatation, if the extensions are infinitesimal, as assumed 
(‘‘very small quantities’). 
2. So far Maxwell. Introducing 


ay 


ey = eg sk ey = e, 
we now consider the case 
en 0.2, = —e,, .€, 0 (2.2) 
and accordingly from (1.13) 
Ps = p(1 mae 2e,) 


Poy = p(1 + 2e,) (2.3) 
Pe ay, 


ll 
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with 2 = E = 
Pn = (Baa + Pos + B..)/3 =o) (2.4) 


This is a case of pure shear.* 

From the Mechanics of Continua we know that if a rectangular coordinate system 
x, y, zis chosen so that the direction x bisects the right angle between the directions 
a and b, while the direction z coincides with c—then the pressure tensor can alter- 
natively be expressed as follows 


1/0 0 O77 0 
By =P |ay—2s 0-170 =p {ay—2s 10 0 | (2.5) 
OF 905-05 62 OF 0 O74) x92 


At the same time we can define a tensor of deformation e,; so that 


C1. <0 Oey ee 
ei; = 0 Co 0 => a 0 0 (2.6) 
Ores Olan g Oe Oe Ones 
This tensor can describe simple shear of the kind 
u,=yy, uy=0, u,=0 (2.7) 


in which 
in whic nee 


is a ‘‘small quantity’’ of the same order of magnitude as e,. Eqs. (2.5) and (2.8) 


then yield - 
f | Bir a) ONY (2.9) 


3. If we consider a deformation of the ‘“‘rectangular vessel’? in which we alter 
the right angle between two sides by a small quantity y, the vessel becomes a 
parallelepiped. If before deformation the edges were of unit length, i.e. the vessel 
was a cube, four of the twelve edges now are of length J1 + 7, which can be 
approximated by 1 + 7/2. 

Now consider a deformation in its reverse sense, i.e. one which deforms the pa- 
rallepiped into the cube, as shown in Figure 1. Let the cube be oriented in the di- 
rections x,y,z, as defined in the preceding Section. We suppose, as before, “‘the 
pressures at first equal (=p) in the three directions perpendicular to the sides’’, 
and enquire after the pressures p,,, Py, Pz, after the deformation. The volume 
of the vessel has not been changed by the deformation, and we can apply 


* Note that in the general case of Eq. (1.13) 
Pm = i, + Pop + p..)/3 = pi = (5/3) ey] (2.4’) 


and p,, is equal to p only if there is no change of volume. 
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Figure 1 


Maxwell’s equation (1.1). To a first approximation e, = e, = e, = 0, as in the linear 
theory, which makes 


e, = 0 (3.1) 


In order to find second approximation terms, we compare Figure 2. The dimensions 
of the ‘‘vessel’’ which have been altered must be taken in directions nor- 
mal to the sides, i.e. in the directions in which the pressures p and p, act upon 
them. We see that only the dimension in the direction in which the pressure p be- 
comes fp, is changed: while the dimension in the direction y, the direction common 
to p and f,,, is not changed. We then find from Figure 1, that in the direction x* 


agat=(1-B) <4 69 


_ * The complete strain-tensor, including the second order term, which has to replace Eq. (2.6) 
is now 


y/2 x2 0 
ey =| 920 : (2.6’) 
0 0 0 


But e, is still = 0, because taking into account the second order term 


e, =(1+ 21 + 41 + 811I)* — 1 


where IJ, II, III are the three invariants, which in this case are 


2 
2g fil bait = 
l= Suto qe HL =0 
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which makes from (1.1) 


ase eat ges (3.3) 


The complete stress tensor is therefore 


Z 


yy 0 
Pij = DP Oy y O 0 (3.4) 
Oo OY 


This implies a second order tension in the direction x, superimposed upon an iso- 
tropic pressure p. I have named such second order stress ‘‘cross-stress’’. In 
general, there will also be cross-stresses in the directions y and z. 

If the shearing tractions t = — p,, acting upon the sides are known, we can intro- 
duce in accordance with Eq. (2.9) 


(3.5) 


AL e 
P 
and therefore from Eq. (3.3) 
2 
% 
une el ae ae 3.6 
ae (3.6) 
4. Now consider the laminar flow of a gas between two practically infinite parallel 
platens at a distance dfrom each other. Let one platen be stationary and the other 
moving with the velocity V. The movements of the molecules will then be composed 
of a laminar motion with velocity 


saad WI (4.1) 


Figure 2 
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Let the line 2~2’-4-4’ of Figure 2 indicate the moving platen, the line 1-3 the 


stationary platen. Then, roughly speaking, all molecules contained at one time > 


between the platens and the planes 1-2 and 3-4 will ‘‘shortly afterwards’’ be con- 
tained between the platens and planes 1-2’ and 3-4’. It may therefore be assumed 
that their thermal motions will cause a stress in accordance with the previous con- 
siderations equal to 


zip tt 0 
Pi; = po;; ed 4 0 0 (4.2) | 
0 0 0 | 


In order to maintain laminar flow in an ambient atmosphere of pressure po, it will 
therefore be necessary to apply to the mass of the gas, in addition to the shearing | 
stress t, transmitted through the two platens, also a tension t”/p somewhere at | 


the free ends of the apparatus in the direction x. This ‘‘cross-stress’’ acting in the 
prism 1-2’-3-4’ will relax and disappear ‘‘shortly after’’ the plane 1-2’ has moved 
into a position 1-2’’, where the length 2’—2” is the same as the length 2-2’, i.e. = yd. 
At that moment, however, another plane 1-2 has moved into position 1-2’, and the 
tensile cross-stress of Eq (4.2) is permanently maintained as long as the laminar 
movement lasts. If the tension t”/p is not applied, the constancy of volume cannot 
be maintained; the mass of gas will contract in the x-direction until p(1 — 17/p”) 
is equal to po, or, assuming po > T, approximately 


2 
a po( 1 ete ) (4.3) 
0 
and 


(4.4) 


5 Now let the two platens be of circular shape with radius R The directions 
r,0,z now take the place of z,x,y, or the stresses ,,, Deg, Pzzs Do, the place of 7,,, 


Pxxs Pyys Pry. The tangential stress j,, = —t is now a function of r, and we can write 
Poz = = tao! Pr = Pyg = 9 (5.1) 
so that the torque is 
R 

M =2n [ cr7dr-= nR*t/2 (5.2) 

10) 

The stress tensor is analogous to (4.2) 

0 0 0 


0 rip rto 
0 IT 0 


Diy = poi; aie 


(5.3) 


| 
| 
| 
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We have to apply the stress-equations in cylindrical coordinates. We assume axial 
symmetry and neglect body forces and inertia terms. The equations are then reduced 
to 


Be 4 Pr sabe Pe | (5.4) 
or 
2 + rt = 0 (5.5) 
with 
PHP = 0 (5.6) 
Integration of Eq. (5.5) yields 
p? = —r’t2 + C (5.7) 


Therefore 
In= p= [C— 178 (5.8) 


As before, we must assume that the gas will contract to some radius rp < R, until 
B,, is equal to po, the ambient atmospheric pressure. Therefore 


C= pp + roto (5.9) 
and 
p= p+ (8-1) <8 a0 
This makes eto -2 
Fo =Ba = ph + (B4rH (5.11) 
and 
Bre OD an —2r % (5.12) 


On Ors ah pees (12! re 


We thus found that when a torque is applied, producing a torsional stress t = rto 
in accordance with Eqs. (5.1) and (5.2), this brings withit a “‘centripetal’’ pressure 
B,,, and an equal axial perssure j,, in accordance with Eq. (5.11). If there is an 
opening at the centre of the stationary platen (r = 0), the apparatus will act asa 
‘‘centripetal pump’’, due to 


D Z D = 4 
Prr| 20 ral = Pool, oo = (D5 ar Foia) (5.13) 
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which is greater than po. At the same time the apparatus forms a self-acting air 
thrust bearing with a bearing capacity of 


To 3 2 3/2 12 
P=2n [ (dre — Po)rdr = a (1 +r 2) —2-3re | (5.14) 
(0) 


0 Po 


These two phenomena, namely the centripetal pumping action and the thrust bearing 
action have been observed in the instruments mentioned above. 

6. In the instruments mentioned, the torque M is not observed, and Eq. (5.2) 
can therefore not be applied. The quantity actually observed is the rotational velocity 
Q produced by the torque M. If the viscosity of the gas is known, we have 


Q 
To = lr (6.1) 


and Eqs. (5.13) and (5.14) can be written in the form 


and 


n°(S) Pin =; E + ron (=) | 32 — Polo HEP (6.3) 


The following qualitative features result from our exploration of the possible 
consequences of Maxwell’s theorem: — 


(i) In accordance with Eq. (4.2), the second order stress superimposed upon the 
first order stress of classical hydrodynamics which appears in plane laminar shear- 
ing flow, is a tension in the direction of flow. In other words, the pressure in this 
direction is reduced 


(ii) In torsional flow, this implies a centripetal pumping effect in accordance 
with Eq. (6.2). 

(iii) The pressure P as defined by Eq. (6.3) is always positive. This implies a self 
supporting thrust bearing action. 

Thesc results agree qualitatively with the observations of Reiner mentioned above. 


Vol. 10C, 1961 M. REINER 11 


7. The foregoing theory is incomplete in the following respect: Maxwell assumes 
that the wall of the vessel is shifted parallel to itself. This would require in our 
case a deformation as shown here in Figure 3. 

We may modify the “‘original’’ shape to the one shown in Figure 4, if we imagine 
every “‘vessel’’ to represent a strip of very small width (comparable with the mean 
free path) and great length d, so that the skew position of the end in the direction 


AY 


cc----- = 
| 
i) 
o 
~ 
+ 
eee original shape we oe ee = OFA! shape 
eee deformed Shope em ees «ere modified original shape 
Figure 3 Figure 4 


y is of no influence. However, we have, in addition, to rotate this strip until it has 
the positions shown in Figure 2, and of this rotation we have not taken account. 
That such rotation materially affects the pressure in the gas was shown by Chapman 


and Cowling. 


8. I am fully aware of the weaknesses in my foregoing argument. I have made 
the present attempt in despair of what I have found in the existing kinetic theories. 
With its weaknesses, my essay at an alternative kinetic theory gives at least the cor- 
rect sign. 

Of objections which may be raised against my theory I mention what Professor 
Cowling has written to me, as follows: — 
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“In your discussion with me you drew the following diagram, 


- ) a+ da | 


and deduced that da/a = 4)”. However, I cannot see why one should not modify 
the diagram to one of the following: 


Figure 6 


The first would give da/a = —4y?, the second da/a =0. These facts convince me 
that Maxwell’s argument is not really applicable to your problem.”’ 

To this I replied: 

“T assume that the shearing flow is from left to right. The fluid approaches in 
a steady state, and the question is after the pressure component p,.,. in the direction 
x, which operates when the sides of the parallelogram are in vertical position. After 
passing this position the state reverts to the steady state with isotropic pressure Po: 
This pressure does not operate in the direction x, but in one inclined to it. In the 
three diagrams shown by you, we have first 6a,/a = }y’, in the second we have 
6a,/a = —4y’, cancelling the first, so that ultimately 6a3/a = da,/a + da,/a=0”. 
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Figure 7 
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A STUDY OF STABILITY AT SUPERCRITICAL SPEEDS* 


HAIM WOLSKI 


Division of Mechanics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


The stability of a shaft rotating at a supercritical speed is discussed. It is proved that two 
degrees of freedom are insufficient for obtaining stability. The effect of friction and the 
behaviour of the shaft in torsion are apparently essential in explaining the phenomenon. 


It is well known that every industrially manufactured rotating part has a certain 
eccentricity causing it to vibrate during rotation. It is also known that at a certain 
speed, called ‘‘critical’’, this vibration reaches maximum amplitude and endangers 
the part by inducing considerable flexural stresses. Above this speed, however, 
motion is steady. A new approach to this phenomenon is presented in this paper. 

Consider a shaft on two bearings, whose moment of inertia is negligible compared 
with that of the central disc (Figure 1). The centre of gravity of the system lies in a 
plane perpendicular to the axis of the bearings as well as to the tangent to the elastic 
curve of the shaft By this means gyroscopic effects are obviated, i.e. the plane re- 
tains its orientation during rotation. 


Bearing 
(Zz) 


FIG. 1 
Description of a flexible shaft 


er This paper was submitted in 1960 to the Technion-Israel Institute of Technology in compliance 
with the requirements for the degree of Master of Science. 
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It is also assumed that the elastic curve is unaffected by the dead weight. Finally 
let us assume an eccentricity e due to faulty production and to non-uniformity of 
the material. The mass of the driving rotor is included in the system, but in order 
to eliminate an additional degree of freedom it is considered infinite. This system 
will be referred to as an elastic shaft. 

It is known that assuming a constant rotational speed @ and solving the equations 
of motion, we obtain: 


—$<—<——<——— 4 
n =A, eo “sin (Jor — t+ ag) 6 a Oe (6t — y) 
Mer = OF + 402 
We nts 2 (1) 
E=A,; e5 "sin (y/o? —wt+ tt) fA eT eae (0t—y) 
/(w2 — 62)? + 4262 
is ‘ 2y0 
y =are 5 2-62 | 


Analysis of these equations shows: 


(1) Friction destroys the collinearity of the points O, S, G, (Figure 2) and an 
angle y is formed. 

(2) After precession at the natural bending frequency has been damped out, 
the centre of gravity describes a circle about the axis of rotation O at a speed 0. 


(3) The radius of this circle is finite for all speeds. 


(4) The phase angle varies, especially about the resonance level, and at a moderate 
supercritical speed this angleis very close to 180°. 


(5) Friction only slightly affects the vibration amplitude, except when 6 approaches 
the resonance level, and mainly affects the lag angle y. 


Schematic diagram of a cross-section through the cg. 
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These conclusions permit definition of the state of equilibrium at a supercritical 
speed (assuming alow damping factor) as follows: 


(a) 9=const. (b) y= 180° (c)q= os e (2) 


The equations of motion, however, fail to provide any information as toits sta- 
bility at supercritical speeds. This aspect was the object of the present study, based 
on the Lagrange equations. 

The procedure is as follows: a slight deviation from equilibrium as defined in 
eq. (2) is assumed; if vibration is found to damp out with time then the equilibrium 
is stable; if it is divergent then, unstable. Let us consider the elastic shaft as a mecha- 
nism with three degrees of freedom (four, not allowing the qualification with regard 
to the mass of the driving rotor). In these circumstances three generalized coordinates 
coordinates suffice for determining the system: the deflection of the shaft (s), the 
phase angle of the centre of gravity (y), and the angle of rotation. 

For convenience, it is preferable to use a rotating system of coordinates with 
origin at S. The disturbed rotation speed 6 + @ should be transferred from 0 to S. 
This imparts to the system a translational velocity (9+ ¢) (s+ 6s) along axis 2 
(see Figure 3). The disturbance ds similarly imparts a translational velocity, s, along 
axis 1. 


(8+ #)(s + bs) 


FIG. 3 
Degrees of freedom ot the flexible shaft 
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The Lagrange formula states that 


d dE, 


= 
Deas ye (3) 


The kinetic energy of a system in combined translational and rotational motion 
is given by 


E, =40 h+4mi?+ ma(p x o,) (4) 


Substituting the following expressions in (4) 


o= O+¢ 
h = (6+ 4)J 
é=p P; = —ecosp Pp, = esinB 
vy =eBsinB+s  v,=eBcosB+(6+ $)(s + ds) 
we obtain 


E, = 4(0 + o)2J +4m {(cisin +s)? + [eB cos B + (0+ )(s + ds)]2 — 
2(6 + ¢)[e2B — e(8 + $)(s + 5s) cos B + essin B] (5) 
For the degree of freedom s we find: 


d OE, 


Ear ots + eB sin B — 0 + )sinf] = m[s + eB(B — 0 — $)cosB + 


e(B — $)sin B] 
Approximating sinf= B, cos B ~ 1 and neglecting all small products such as 
| B2, Bo, Bo etc., linearisation is obtained and we have 


d OE, 


i 3 = m(s - ebp) (6) 


or, substituting q = s—e 
Es — m[(o + f(s + 83) + ef(0 + G)008 f - €(0 + 4)# 008] ~ (7) 
m(q@2 + 2g0b + 626s + e0B) 
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At this stage we consider the following three factors: 
1) Flexural elastic recovery force P, =cs 
2) Friction, proportional to the rotational speed of the c.g. P2 = UmOG 


3) Torsional elastic recovery moment M,=k¢ 


There are additional friction forces counteracting any change of coordinates. | 
These will not be allowed for at this stage, although any withdrawal of energy from 
the system favours stability. We have 


OX, _ z x5 d., 0x3 
eee 
Ox 
P, ae —c(s + ds) (8) | 


Combining eqs. (6), (7) and (8) we obtain the first equation of motion 


m(s — 2e08 — 40 — 246 — 675s) = —c(s + 4s) (9) 
or 


5S — 075s + w25s — 2e08 — 20qd = —sw? + gO? =0 (10) 
A Laplace transformation, denoted by a bar, will give 


(p* — Q”)5 — 2eOpB — 20qp¢$ = sp (11) 


For the generalised coordinate y we have 


ay dae [(B — dye + sh + 05] (12) 
a ~ meO(sOB — Ss) (13) 
Ox, OX : 0x3 ee 
ag ape aca 
0 
M, aE = —1,,eq0 (14) 


Combining eqs. (12), (13) and (14) we obtain the second equation of motion: 


eB + qb + 20s — s02B = —2ng0 (15) 
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| 


) which, transformed, gives 


20ps + (ep” — s6?)B + qp?B = —2ugO/p (16) 
For the rotational degree of freedom, 


d dE, 


“Oe aaa ca, (17) 
OE, 

ape (18) 
OR Pe ONE nat | Oe en 

Ogg sch gaa dota 

>, = He g(O + ¢) — ko (19) 


Combining eqs (17), (18), (19) according to the Lagrange formula, we obtain 
Jb + mq(20s + eB) = —kb — ppeg(O + ¢) (20) 
After the following substitutions k/M = w?, J/M = i? 
i2g + 2q0s + eqB +W2o+ 2ueqd = —2yeq0 (21) 
transformation gives 


2q0ps + eqpB + (i?p? + 2ueqp + W?)b = —2yeqg6/p (22) 


We now havea system of three equations in three unknowns, eqs. (11), (16), (22). 


_ The roots of the principal determinant of the three equations determine the stability 
_ of the system. If their signs are negative the system is stable. 


The principal determinant is : 


p?>-Q? = —2eOp —2q6p 
6 
20p ep?—s0” gp” = 2 a;p 


e 23 
2q0p eqp’ i?p? + 2ueqp +" ce 
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Still another substitution is made in the polynomial obtained, namely 


K = (3- =) +0; (24) | 


and assuming further as in engineering practice that i? > q* we determine the co- | 
efficients of the polynomial as follows: 


2 


Ag = ei 
a, = 2ye*q 

a4 = ew? +i7eK 

a; = 2pe*qK (25) 


a, = weK + ei0* 
a, = 2ue?q0* 
ey70* 


ao 


By the Hurwitz criterion: 

1) The coefficients of the polynomial should be positive. As for the terms con- 
taining K, we see that for 0 0, se, hence K >0. The transitional value, for 
which K = 0, should now be found. Then by (24) 


(3 — s/e)@? +a, =0 | (26) 


2 


Beri Wee Sea tee oe } 4 ; 
substituting > = Gn Bives the biquadratic equation 


20* — 2w20?— wt =0 (27) 
The only solution of physical significance is 


0=1.170, (28) 


This result is quite logical since the calculation is based on the assumption of 
y = 180°, which is valid only above a fairly high supercritical speed level. 

Friction also has a certain influence, but as it had been disregarded in eq. [2(c)], 
it does not appear in the result. 
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2) The following determinant, and all minors, should be greater than zero. 

as a3 : ay 0 0 0 

a6 a4 : a2 | ao : 0 2 0 

0 as a3 | ay | 0 | 0 

0 Ag a4 Gee U code on 0 
) 0 0 as a3 ay 0 
| 0 0 a6 a4 a, ao 
| Condition (a) 
) asa, — aga, > 0 (30) 
_ substitution of coefficients gives 
2ue*qy? > 0 (31) 
Condition (b) 
) A3(45a4 — 463) + as(aga, — asa) > 0 (32) 
_ Here substitution of coefficients actually gives zero. 
- Condition (c) 

(4544 — 643)(a342 — 44a, + Ads) — (aya6 — 45)" >0 (33) 


Here again the expression vanishes on substitution of coefficients. 
Condition (d) 
(42 4, — 4349) [(a5(44a3 — a5a2) + a6(2asa, — a3?)] — a5(aoa5 — a4a4)?— 
— 46434;(4o45 — 4444) —agaj > 0 (34) 
Same as above. 


Condition (e) is obtained by multiplying eq. (34) by ao again giving zero. 


CONCLUSIONS 


1) On the above assumption, expansion of several successive minors gave zero, 
indicating that the roots are on the imaginary axis, but owing to linearisation their 
exact location cannot be ascertained. In practice this calls for iteration with a higher 
degree of accuracy, but it can be assumed that in view of additional damping (for 
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example, in each degree of freedom) the system is certain to be stable. It is clear, 


in any event, that an elastic shaft is a borderline case from the point of view of | 


stability. 


2) It appears that friction on one hand, and behaviour in torsion on the other, 
are of decisive importance in this respect. The elasticity of the shaft in torsion intro- 
duces an additional degree of freedom. It will be proved later that two degrees of | 


freedom are insufficient for obtaining stability. 


3) When friction is not high, 0 = 1.17 @, is the transition point between stability | 
and instability for a 180° lag. For lower supercritical speeds, smal'er angles should | 


be taken into account, namely 90° < y < 180° which tends to complicate the cal- 
culation. 


4) It is unnecessary to postulate a linear relationship between friction and speed. | 


If square-proportionality is assumed, eq. (19) is modified as follows: 


Ox 2 é 
—M, ap = Um (8 +)? & lmqeO? + mg Bed 


The determinant (23) will include an element containing o differing from its pre- 
decessor by the factor 20, which does not affect the results obtained so far. 


5) As already mentioned, it can be shown that two degrees of freedom are insuf- 
ficient for obtaining stability. In order to verify this, the order of the determinant 
(23) should be reduced by one by eliminating the rotational degree of freedom. 
The corresponding minor is 


p? —Q? —2e0p 


= e(p* + Kp’ + 6*) (35) 
20p ep” — s6” 


Stability is conditional on p being imaginary and p” real and negative, 


2p? = —~K+t x — 46? (36) 
as well as on 
a) K>0 satisfied 
b) 467>0 satisfied 


(37) 


c) K*>46? not satisfied for any 
supercritical speed 
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In these circumstances, even if friction is reintroduced for each degree of freedom, 
stability cannot be obtained at any speed. This refutes the hypotheses of Robertson’, 
Den Hartog? and Stodola* on stability under two degrees of freedom. 


6) If k in eq. (25) is allowed to tend to infinity, the following expression is ob- 
tained: 
ye(p* + Kp” + 0*) (38) 
By (37), it is seen that a torsionally rigid shaft is unstable. 
For k =0 no conclusions can be drawn owing to linearisation, but it seems that 
there exists a transitional value of k between stability and instability. 


NOTATION 
ce = flexural ‘‘spring constant’’ of shaft v,,v. = components of v, along axes 1 
E, = kinetic energy. and 2 resp. 
e = eccentricity B = variation of the phase angle 
G = centre of gravity y = phase angle 
h = moment of momentum 6s = changes in deflection of shaft 
J = moment of inertia of disc about 4 = abscissaofcentreof gravity with 
the geometric centre. respect to fixed system 
i = radius of gyration (i= ./J/m) 
k = torsional “‘spring constant’ of shaft @ = rotational speed of shaft 
m = mass of disc i sem 
M, = moment of internal force Lm = damping factor 
O = centre of rotation € = ordinate of centre of gravity 
P, = internal force with respect to fixed system 
p = Laplace transformation p = position vector from origin of 
q = radius of centre of gravity variable system to centre of 
q; = generalised coordinate gravity 
>= centre of disc P1,P2= components along axes1& 2 resp. 
s = deflection of shaft go = increment in angle of rotation 
; -= time due to speed changes 
v, = translational velocity of centre of @ = rotaticnal speed in general 
gravity ©, = natural bending frequency 
Oe 
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ENERGY LEVELS IN MODES OF FREE EXTENSIONAL VIBRATIONS 
OF A FINITE PLATE* 


AMNON Foux AND NoRMAN DAVIDS 
Division of Mechanics, Technion-Israel Institute of Technology, Haifa 
and 


Department of Engineering Mechanics, The Pennsylvania State University 


ABSTRACT 


Energy requirements for exciting free extensional vibrations of finite plates are analysed 
in this paper. The solution for an infinite plate is utilized and a superposition method 
applied to account for the free periphery of the finite disk. An illustrative numerical example 
presents the relative magnitudes of the required energy for the first three modes of 
vibration. It is found that the resulting energy ratio for the third mode is less than that for 
the second mode. 


NOTATION 
a = one half the thickness of a plate v = volume 
A,B,D = arbitrary constants W = strain energy density 
ty. g ae 
h Soe ei aes pyr = W,. work done ona cylindrical 
A+ 2u face in an infinite plate dur- 
J, = Bessel function of order n ing one fourth of a period 
2 2 
k =a parameter = /02 2 ) = wave number = = 
M, = functions of y,h,k and a A = dilatation 
p = frequency x 2x & = strain 
r,0,z = cylindrical coordinates A = wave length 
R, = radius of finite plate 2,4 = Lame’s constants 
5 = surface v = Poisson’s ratio 
t = time o = stress 
u; = displacement along the ith axis t = period 


U,U’ = strain energy 


* This work has been partially sponsored by the Office of Ordnance Research (now A.R.O.D.) 
U.S. Army and The Pennsylvania State University. 
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1. INTRODUCTION 


This paper presents an energy analysis of flat plates subjected to free radial vib- 
rations. The investigation was prompted by observing the situation for buckled 
plates, in which some higher buckling modes may occur having lower energy states. 

It is generally accepted to associate the frequency of different modes of vibra- 
tions with thei: wave length, as in the case of axial vibrations of rods, where higher 
modes have higher frequencies and shorter wave lengths. In the case of vibrations 
of plates where th2 displacements are not plane, modes are possible for each major 
mode associated with a certain wave length. The difference in these modes would 
be in the configuration of the displacements as a function of z and not in the wave 
length. In the present problem, such a comparison of energies for the different modes 
would afford some insight into certain questions of stress-wave propagation in flat 
plates. One such question concerns the relationship between the type of excitation, 
and the mode, or combination of modes, which is actually excited. Another type 
of question concerns the interpretation of observations, such as velocities of pro- 
pagation, or natural frequencies, in terms of modal behavior. 

In both these problems, the relative energy levels among the different modes 
play a significant role. In the present paper, although the problem is not considered 
in full generality, an energy analysis is made for the finite plate for the first three 
modes. It is based on the solution obtained from the infinite plate excited at the 
center. The legitimacy of this procedure, even though a mathematical proof is not 
given, is suggested by the fairly good agreement between the analytical solution 
for infinite plates given by Davids and Kumar! and the experimental data of fiee 
vibrations of discs from the Bell Telephone Laboratories given by Kane and Mindlin’. 

It is found, for the case where the radial amplitude at the central plane is kept 
the same for all modes, that the energy required to excite the third mode is con- 
siderably less than that for the second mode, although both are higher than the 
first mode. This would mean that the third mode is here easier to excite than the 
second. The pulse shape, however, is actually the determining factor. 

In Mrs. Feder’s3 experiments at Frankford Arsenal it was found that the velocity 
of stress-waves was substantially higher than that calculated directly from the elastic 
constants of the material based on the lowest mode. It is not known at present 
what the actual reason for this is, or whether the results of this paper are directly 
applicable to the conditions in her experiments. Nevertheless our results suggest 
the possibility that higher modes are being excited, for if so, higher observed velo- 
cities would be interpreted differently. 


2. EXTENSIONAL VIBRATIONS OF A FINITE PLATE 


The simplest form (lowest major mode) of free radial vibrations of a finite disk is 
the case where it has a node at its center and the next antinode is at thz free pevi- 
phery. The total energy in the disk, which is constant, is the sum of the kinetic and 
the strain energy, and these are being transformed from one state to the other. 
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For the case of harmonic displacements with time, the maximum strain energy 
is equal to the maximum kinetic energy and these are equal to the energy originally 
transmitted to the disk in order to bring it to its state of vibrations. 

Free vibrations of infinite plates are incomprehensible in the simple sense of 
modes of vibrations, since they are functions of both the frequency and the wave- 
length. This property enables us to choose the wavelength such that the portion 
of the infinite plate of radius R, = A/4 simulates the finite disk (Figure 1) with the 


SPL LF LE LP LLP LP PAP LPP LP LPP 


2R, 


a a II IPF 08) EE ES 


Figure 1 
Correlation between the finite and infinite plates 


exception of the free periphery possessed by the finite disk. This however can be 
accounted for by considering the work done by the restraining stresses on the 
radial boundary of the cut section of the infinite plate. That is, the maximum strain 
energy in the finite disk is the sum of the maximum strain energy in a portion of an 
infinite plate of the same radius, and the work done on the radial boundary of this 
section by the internal stresses (Figure 2) during the time t/4 whichis the time required 
for transforming the kinetic energy into strain energy. 


Figure 2 
Stresses on a radial boundary of an infinite plate 
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3. STRAIN ENERGY IN INFINITE PLATES 


The strain energy U is defined4 as 
U = f Wdv (1) 
where Wis the strain energy density, given in tensor notation ,by 
W = to; ,8;; = 4AA? + pez, =) 
Thus for the strain energy we have 
eee fidtde . f 2ue?, dv (3) 


where A is the dilatation, A = (¢,, + &@ + €,,) and A and p are the Lamé Constants 
of the material. 


When equation (3) is written in an unabridged form it appears as 


Uw 44 [A2dv + 2u f (2 Tv Ge pebie 2 eoae, FO Dee 2¢?,)do| (4) 


The displacements are given! as: 


A 


u, = (Acoshz + Bcoskz)J,(yr)e" 


wey Us B (5) 
u, = (=-sin hz —- —* sin kz ) JoCvre'”* 
The strains are defined* as follows: 
Ou, 
Epp = or 
1 Ou, , &  u, 
9 SeMA0y | aes 
pe 
ee) 


AiR pie 
aes) 
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Since 3 1 3 
ap LIiOr)] = Wolvr) — iQ); 3 Joly) = — sr"), 


We have for the strains 


(Acos hz + Bcoskz)[yJo(yr) — a yr) ]e?" 


Epp = 
1 ipt 
&99 = (Acoshz + Bcoskz) . J ,(yr)e” 
2 
&, = (= cos hz — By cos kz) Jo(yrjei”" (6) 
Eg = & = 0 
= é B 7 F ipt 
&, = [- Ahsinhz + =(5- k) sin kz | J, (yrje 


For the dilatation we have 


A =\A (‘ = r) cos hz Jo(yr)e”* (7) 


The terms of equation (4), which represent the strain energy in the plate at any time, 
are complex and appear as squares. Although complex quantities are employed in 
the analytical operations, the physical entities must be finally represented by real 
functions. In the case of expressions such as (4) involving squares or products, we 
must multiply the real part of the factors, instead of taking products. 

The squares of the real parts of the complex quantities are these of the coeffi- 
cients of e”’. This leads to the expression for maximum strain energy in a section 
of radius R, of an infinite plate as follows: 


a 


UL = ; [2 af i) (4 (< — 9) coshe Jo(yn} rdr dO dz (8) 


Caan’ § 


+ auf f ff [[(Acos he + Beoske) (rar) - 41009)} 


fe {(4 coshz + Bcos kz )* J,(yr) \'4 + (A eos hz— By cos kz)Jo0n)} 


+ 2{(— Ahsinhz + ze — k) sinkz ) AQ) | rdr dO dz 


where U’ denotes the strain energy in a finite portion of the infinite plate. 
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Performing the integration of (8) with the substitution of 


_  2Ahk sinha 
~ (v2 =k?) sinka (8a) 
as obtained from equation 14 of Reference 1, and 
IRs 5 
TO oe 
as explained above, the following expression for the maximum strain is obtained. 
} Ara? 1 1 1 2 
Oak y2 [Mig Dm | + o[Mar? ( ILE taal ea rary 2 ms) 
1 1 
+ M35) m, + Mas) ma] } (9) 
where 
es (ae! OC 
ym ae? (ms (me iyi = 04439 
S (=1)"2m + 1)1(n/4y"*? 
mM, = SN = 0.4749 
yma = de mile FIP (10) 
2 (=1)"(2m + 2)!(0/4)?"*? 
= EES SE Nosh ed = 0.4551 
Doms Mto [(m+1)!]2(m +2)! 
-y (-1)"(Qm+ Gia 0.2028 
dm = m=o Mm! (m+ 1)! [(m 4 2)!]2 : 
and 
h?—y?), sin 2ha 
Biel) (42 a) 
sin 2ha 4hk sinha jsin(h—k)a_— sin(h+k)a 
My =} (04 2h beer sinka ( h—k— = hk 
2hk  sinhay? sin 2ka 
ua & — k2 aia (4 2k )} 
h* sin 2ha 4h®k sinha jsin(h—k)a_ sin(h+k)a 
= —? nee eS ee a ee Ss - a SS 
Pica anergy ca ae 


(i ta (oe 


me sin2ha, __, sinha jsin(h—k)a_ sin(h+k)a 
My = iY (a ~ on) Anka | (h—k) (h +k) 


eis) Cena 
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Thus we finally have 
A2 


Umax = 


3 
a [A(0.0555 M,) + 1(0.0520 72M, + 0.1110M, + 0.1014 M,)] (12) 


4. WORK DONE ON THE RADIAL BOUNDARY OF RADIUS R, OF AN INFINITE PLATE 

We now consider the effect of the work done on a cylindrical boundary of radius 
R, of an infinite plate. Since this is a contribution from the radial stresses and the 
vertical shear stresses during a quarter of a cycle, the expression for the work is then: 


t/4 1/4 : 
We= f fo,ujdtds= f f (ou, + o,.U,)dtds 
OF us 0 « 
or 
oa . . 
W, = f ff [(AA + 2u8,,)u, + 2u8,,u,] dt ds (13) 
ORs 


It is noted from equations (5) (6) (7) that the terms inside the integral are products 
of complex quantities of the form C, e'”’ x ip C, e'”' where C, and C, are real func- 
tions of the coordinates. 


On integrating with respect to time we must consider the product of the real 
parts only. 


The real part is given by 
7/4 
f Re(Cye )x Re(ipC,e)dt- 
(0) 


1/4 
= 4 f (ipCyCre * ipC,Cre ~** ipCyCz + ipC,C2 )dt= Se: (14) 


The work then is one half the surface integral of the product of the coefficients of 
e’' in the strains and the displacements. 
Then, 


2n a 


W, = asf [ {a4 (ea cos hz Jo(yr) + 2u(Acoshz + Bcoskz) 


(»Jo(or) - *J,07))| (Acoshz + Bcoskz) J, (yr) 


+ 2u (- Ahsinhz + ie : - sin kz) J,(yr) 


Ah ., By , 
(Se sinhz — =, sin kz ) Jo(yr) | rd0dz (15) 
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Integrations with respect to 6 and z for constant r= R, =(z/2y) and substitution 

2hk sinha 
(p= A 

. y2 — k?2 sinka 

ary per one fourth of a cycle: 


w= A [a bro(S):(3)} me: 


yields the following for the work done on the radial bound- 


) 
Hl? ErelF)a(5) - 8G) +270 (F)(F)me}] 8) 
: where 
Po fae : ; Se pee : 

p= 2 es St) eae 

Zz * 2 Z 2 * 
= et ee Oe Be m 
| 
pein Hea pial ale) oc 2 aay IDE) SI ey 

s +k y2 — k2 \sinka 2k 
and M, is as given in equations (11). 
Finally we have 
W, = ae [A(0.2675 Ms) + (0.1259 yM, + 0.5349 M,)] (18) 


5. MAXIMUM STRAIN ENERGY IN A FINITE DISK 
The maximum strain energy in a finite disk is the maximum strain energy in the 


portion of an infinite plate of the same radius plus the work done on the periphery 
of this section by the internal restraining stresses during a fourth of a cycle. That is 
Unax = Us = W,, 
Thus, combining the expressions (12) and (18), 


Vee ine 7 [40.0555 M,) + u(0.052072M, + 0.1110 M; + 0.1014M,)| 


oe [4(0.2675 Ms) + 1(0.1259 yM, + 0.5349 M.)|} (19) 


where the M’s are given in equations (11) and (17). 
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6. CALCULATIONS FOR THE ENERGY-LEVELS OF A STEEL DISK IN ITS FIRST THREE MODES 


The strain energy in a plate is a function of the mode of deformation of the material 
and the amplitude of the displacements. For the purpose of comparing the maximum 
strain energy in the different modes of vibrations, the amplitude of deformation 
must be kept constant for all modes. It would be reasonable to fix the radial dis- 
placements at z = 0 on the periphery of the disk for all modes. Thus, as can be 
seen from eq. (5) and (8a) 


2hk ~=sinha 
A+B=A(1 + Gf =) | =D) = COnSt: 
or 
D 
ae 2hk sinha (20) 


1+ 


y2 — h2 sinka 


For an arbitrsrily-chosen steel disk with dimensions a = 0.324cm and R, = 0.6707cm 
we have y = 7/2R, = 2.342. From the velocity curves! the values of p for the dif- 
ferent modes can be obtained. Subsequently the h’s, k’s and the M’s are determined 
as shown in Table I. 


TABLE I. 
y = 2.342 
1st Mode 2nd Mode 3rd Mode 
p2 15.96217 x 1011 72.56386 x 1011 131.28361 x 1011 
h i 1.0375 3.8153 §.5474 
k 3.1536 8.0412 11.0192 
M, 5.2824 6.0806 33.1712 
M2 1.4851 0.7084 3.1419 
M; 1.7592 26.7182 40.9023 
Ms, 0.0304 6.9653 93.4594 
M; —2.7946 0.2429 6.6288 
Mg 0.0132 1.8140 —2.5034 


On substituting the numerical values for the first three modes in the expression 
(19) for Ung, we have for the maximum strain energy 
(Umax)r = 2.2459 x 10'* (nD)? —_ for the first mode 
(Unmax)r1 = 30.2608 x 10'*(xD)* —_ for the second mode 
(Umnax)111= 8.3701 x 10'*(2D)? for the third mode 
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FIRST MODE 


NODAL PLANES 
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SECOND MODE 


THIRD MODE J 


Figure 3 
Radial displacement at plate periphery 


The relative energy needed to excite the second mode with respect to the one 


needed for the first mode is, - 
(Umax) 1 (21) 
See 13.5 
(Omax)r 
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and the one of the third mode to the first mode 


(Unadur ah 
eae a7 (22) 


The configuration of the radial displacement on the periphery of a disk (Figure 3) 
makes physically plausible the result that the energy needed for the second mode 
1s higher than that needed for the third mode. 


7. GENERAL CONCLUSIONS 


We are not attempting here to find the values of the absolute energies in finite plates | 
but only their relative magnitudes. For this reason the use of the infinite plate solution 
may be assumed to be of the same general nature. Although this is not proved mathe- 
matically, it is suggested by the experimental agreement, as mentioned before. 
Furthermore, a calculation was made by taking the wave length as for harmonic 
function of r instead of the Bessel function and it was found that the values of the 
total energy differed by about 20%, but the relative energies are of the order of 
5% difference. It can also be concluded that even though the relative energy level! 
for the different modes, associated with the longest wavelength, was discussed in the 
previous sections, it is the belief of the authors that the pulse characteristics are 
the determining factors for the mode excited. This last statement is limited to the 
case where the pulse shape corresponds to the deformation configuration of a single 

mode, otherwise a complex combination of modes may be generated. | 
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ON THE DYNAMICAL DERIVATIVE OF PHYSICAL TENSORS 


ZVI KARNI 


Division of Mechanics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


The spatial concept of the curvilinear derivative of a physical tensor is extended to 
include also the time derivative in a space-time manifold. 


0. The summation convention adopted here implies that repeated Greek indices 
only are summed whereas repeated Latin indices should not be summed. Cartesian 
coordinates are denoted by y;, curvilinear coordinates by x;. The discussion is con- 
fined to orthogonal coordinates hence all indices are written as subscripts. 


1. Consider a continuous and single valued vector field A with respect toa three- 


dimensional curvilinear orthogonal network formed by the three families of curves 
x; = const . (i = 1,2,3), where 


Vi = Vi(%1, X25 %3) 
The Cartesian coordinates are the components of the position vector r 


(tr), = 


not so with the curvilinear coordinates x;. 


At any point P of the network, a local system of reference is formed by the unit 


| tangents 
ct or 
ag | Os; 
of the coordinate curves where 
ds, = J 9: 4x; (1.1) 
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and g,; are the components of the two-fold metric system (which is not a physical - 
tensor) found from | 


The local system is thus a Cartesian system, and when the network is orthogonal, 
it is orthogonal-cartesian. 


The projection of the principal curvature 1/p; of the coordinate curve r(x,,a,b) 
(a,b—constants) on the plane t, x t, = — t, x t; of the local system at P is defined | 
as the physical Christoffel symbol G;,,. It follows from the definiticn that 


1 


Pijk = Gin = — Ging (1.3) 


We also find that for three mutually different indices 


because the principal curvature 1/p; has no projection onits tangent t;. In view of 
(1.3) and (1.4) the matrix of the non-vanishing symbols assumes the form | 


0 Git Gis, 
[Gui .23= —[Guj] =) Gai2 0 G232 (1.5)} 
G313, G323 0 | 


It can be shown* that the G,;, are found from the components of the metric system 
gi; by means of the relation 


1 Ogun 09; | 
C= & ~ 3) (1.6); 


ay) GiiD iI 


2.The physical components of A at each point P of the network are identified with 
their projections on the axes of the local system at that point. Thus, and if in order: 
to illustrate it better we restrict ourselves to the plane as shown in Figure 1, the) 
components of A are (Figure 1b) 


* A more detailed discussion on the physical Christoffel symbols can be found in a paper by 
Braun, Karni and Reiner to appear shortly. 


a 
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Figure 1 (a) Figure 1 (b) 
A vector field with respect to a curvilinear The curvilinear derivative 
network 


(A), =A, ba PK; (A), =A, a PL, 


those of A + fli (ds, = 1) 
Os, 
Ye eae = 
1 


It is readily observed that the scalar differences of the projections namely 


PM—PK; PN-PL 


; 0A 
are not equal to the components of the vector difference mes because the latter are 
1 


represented (if we neglect differentials of the second order) by 


OA =— OA —— 
SS = = — = AS 
(Fe), aS (3), 


? oA : 
The physical invariantis obviously the vector difference a and we refer to it as the 
j 


‘“‘curvilinear’’ derivative of the vector field. This derivative is denoted by a comma 


thus aA 


A, = (2.1) 


as, 
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From Figure 1b it is observed that for a unit arc length ds; 


= PM ae) 
(5) - KR = PM cos(day21) — PK + PN sin(daya1) = 5° + A2Gi21 
Lal 
also 
a 0A 
(3) = LS = PN cos(da,21) — PL ~ PM sin(diy21) = 35? + ArGur2 
1/2 


Or in general the i-th component of the curvilinear derivative, denoted A; ; is 
given by 
OA 0A, | 
A, ; =(A,,))i = (5) - as + Ae Gig Qy 


This result, although developed for the plane, applies equally well for three (and | 
more) dimensions. On extending further the curvilinear derivative of vectors to 
absolute tensors of rank n by means of successive outer products, one finds 


OAs, t2...in 

iteints = isqhe Pt Agis..cin Gyat, + Ajo .ntaG jat, + c--Aigig + oF jat (2.3) 
The aes derivative is a partial derivative. On multiplying Eq. (2.2) or (2.3) 
by —— <2 and contracting on j, we obtain the total derivative along a curve of a vector 


or a tensor field, which is denoted by the variation symbol 6. Thus 


5A, dsy_ aA, ds, 
Oy a, 7 gee ae G4) 
OAs tn, CAlinate ds 
5s ae oi ete (Ag) 4, Goat, ms Ane, Gpai, a8 tas tenad Goat, a (2.5) 


The total derivative of physical tensors is analogous to the intrinsic derivative 
of the ordinary, analytical tensors. 


3. Referring back to Figure 1b, one realizes that the appearance of the physical 
Christoffel symbols in the expression for the curvilinear derivative is due to the 


~~ —— << ) 
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fact that the local systems of the network are rotated each with respect to the other, 
the projections of the rotation on the planes of the local system being denoted by 
da; ;, which are equal in value to G;;, for unit arcs ds;. Now, it is immaterial 
whether the local systems happen to be rotated or are being rotated deliberately 
by say a rigid rotational motion, pending the angles of rotation are the same. 
Therefore, if instead of the angles per unit arc 


Ae; i, = Win; 


ds; ce Gin (i,j,k = 1,2,3) (3.1) 
the angles per unit time, i.e. the angular velocity components oj, = —@,j;, are in- 
serted into Eq. (2.2), the following expression for the i-th component of the vector 
difference per unit time is obtained 


dA dA, 
(4t)- ap a he A,Mai (3.2) 


Equation (3.2), which may be regarded as the ‘‘dynamical derivative’’ of a vector, 
can be combined with the curvilinear derivative so as to form a single “‘physical 
derivative’’. For this purpose we consider a space-time continuum with four coor- 
dinates x;(i = 1,2,3,4) and let x, be the time coordinate, hence also 


dt = dxX4 = dy4 = ds4 


Next, the components w,, are rewritten as three index angular velocity components 


dt, ; 
Oj, = Ws jn = a (3.3) 


Equation (3.1), consequently also Eqs.(2.2) and (2.3), now hold for the four-dimen- 
sional continuum with the further properties 


Ga jx = O jx (j,k = 1,2,3) 3.4 
G44, = 0 (k= 1...4) oe 


The matrix of the spatial Christoffel symbols — Eq. (1.5)— is now reorganized 
so as to form with relations (3.4) a single matrix for the non-vanishirg physical 
Christoffel symbols in a space-time manifold 


1) 
40 Z. KARNI Bull. Res. Counc. of Israel | 


1 1 
0 G 0 
Giz sty P112 P131 
1 iy 
Gaiz Gr23 0 P212 + P223 
Giixl] = — |Gajl = = ; (3.5) 
: C223 Caa1 P323 P313 
G4i2 Gazz Gast O12 W230 O31 


In this picture the angular velocity components of a rotating local system form | 
part of a pattern of components of a physical tensor of rank 3 in space-time—the | 
physical Christoffel symbols—although in practical applications only three non | 
vanishing components remain which make the angular velocity look like a vector; 
namely it is in fact a pseudo vector. 


4. As already mentioned in Sec. 2, the physical derivative of a vector and a 
tensor is a physical invariant in respect of a transformation from one local system 
to the other. Thus a differential vector equation formulated in Cartesian 
coordinates is automatically expressed in any curvilinear network if each of the | 
rectilinear derivatives e is substituted by the corresponding physical derivative 
A, ; (similarly by the total derivative along a curve (the intrinsic derivative), because 
the two in fact represent the same thing. The same applies to a differential tensor 
equation of any rank. This procedure offers a wide range of applications as the 
transforming of a Cartesian tensor equation to curvilinear orthogonal coordinates 
is now achieved in a routine and straightforward manner. As examples, consider 
first the acceleration a of a particle in space. When referred to a Cartesian 
network, we have 


_ a de 
dt dt? 
or in terms of components 
dv, d?r; 
Oy = as (4.1) 


It should again be emphasized that the position vector of the particle ris a physical 
vector only in a Cartesian network because in this case its projecticns on the re- 
ference system coincide with its components. With respect to curvilinear networks, 
only its first and higher derivatives are physical tensors. Thus, and in accordance 


3 ye 
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with what has been said, the general expressicn for the acceleraticn cf a rarticle 
with respect to a rotating system in a rectilinear (Cartesian) network reads 


6°r,; 67, 
oF = yh pe 
and in curvilinear networks 
_ by _ ov; 
Tie as, Ge) 


In the expansion of Eq. (4.2), the expression for the second curvilinear derivative of 
a vector is needed. This is found from Eqs. (2.2) and (2.3) 


aA, @ A, deg , OA 
Ain, = a + Aj,qGraj + Ax,jGrai = as, Pele oe jai t ao Gra; 
0A, dG ai 
+ Jy Guai + Aa Got Apl Gopi Gras + GipeGrsi)- (4.4) 
With 
0 0 0 
0 0 0 
[|] = lol=(4]3 Ga = Oe enlace 
v, 


Dagnre By in is 


Eq. (4.2) reduces to the evaluation of r;,44. Setting j = k = 4 and r; for A; in Eq. (4.4), 
also denoting by a dot the derivative with respect to time, the following familiar 
expression for the x-component of the acceleraticn for example (i.e. for i = 1) is 
obtained 


a, =x + Aw,v, — w,v,y) + ZH, — yO, — x(0,? + @,? + w,7) + @,(x@,+ yO, +20,) 


the second term being the Coriolis acceleration. 


§ For a rotating system in the plane of a polar-cylindrical network, we have 


O10 0 

r by dr a: , id ob 

[>] = 3 |ds;] = ; jal = ai 5 WG =| ae 
¢ Ge | @.’ 0:20 
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Equation (4.3) now yieldsin view cf Eq. (2.4) 


r —1r0* — ro, 


a, = 
ag = 270 + 70 + fo, 
a,=2Z 


The flow tensor f;; is defined as 


: 1 0 fou; . Ou 
fy = 4; = 5x (Ft os) (j= 123) 


where e,, is the strain tensor and u,(y,,y2,)3,) the displacement. In view of Eqs.(2.2) | 


and (4.4) the following expressicn for the ccmponents of the flow tensor with respect 
to a rotating system in a curvilinear network is obtained 


0 Ou; 6 Ou; Ou; 


2f;; = c. (u;, J4 ae U;, i4) = 03; Ag = 6s; Ct. eos a Bp (Gazi na Ginj) a = as, att as, Gua 
au, du 
1 as; rege as 0s; a Gaui + Ugr, ot vay r G jai) or Ug( Gag.Gaaj a5 Gap ;Gaai T G jpaG gai + GipaG aaj) 


For the polar-cylindrical network mentioned above we find 


_ OU, OU OU pibiivthe 
In = Ge (35+ 5-H )o, 


_ Oy | ty Guy Clg. Ue 
foam Sit Sea Sie 4 it — St) w, 


| ; 
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UTILISATION OF DIGITAL COMPUTERS IN STRUCTURAL ANALYSIS 


J. LOEWINGER 


Building Research Station, Technion-Israel Institute of Technology, Haifa 


I. THE SIGNIFICANCE OF THE COMPUTER IN STRUCTURAL ANALYSIS 


Efficient use of the electronic computer in structural analysis necessitates a new for- 
mulation of the theory. In this respect, the following characteristics, distinguishing 
the computer from the ordinary desk calculator, must be taken into account: 


electronic speed, 

infallibility, 

programme stored in advance, 

capacity for independent modification of instructions. 

capacity for logical decisions (e.g. discontinuation of process, choice of cross- 
section," cic. ), 


These properties permit recourse to exact analytical methods hitherto considered 
impracticable including such complications as variable moments of inertia, non- 
linear stress-strain relationships, nonuniform sections, settlement of supports, 
displacement of joints, and curving of plane sections under deformation (€.g. in 
curved bars etc.). By introducing a variable modulus of elasticity, instability pheno- 
mena in the plastic range can also be analysed. 


II. SPECIAL REQUIREMENTS AS TO THE TYPE OF COMPUTER 


The computer comprises the following components: 


are Control —er Te 
| Input BeaTS= . eo SS 
Output Process 


Individual computers vary mainly with respect to the input-output units and the 
type and size of the memory device. 
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Punch-card input-output units are best suited for the problem in question. 


Formulation of the problem is governed, in some cases, by the available memory ke 
capacity. High-capacity computers are especially convenient for matrix calculations; | 
it can be shown how low memory capacities can be used to satisfaction. The respec- | 


tive advantages and disadvantages of the floating- and fixed decimal point arith- 


metic devices and the available subroutines and interpretation programmes should | 


be considered. 
Ill. THE PROCESS OF STRUCTURAL ANALYSIS 
1. Formulation of problem (suitable basic systems, etc.) 


2. Mathematica’ formulation of solution 


The advantage of the matrix formulation of problems in structural analysis lies | 


in its universal suitability for programming. Apart from this technical advantage, 
matrix formulation permits simple resolution of space and plane systems into suf- 
ficient numbers of linear components (e.g., slabs ahd shells into beam- and arch 
grids resp.) thus reducing their analysis to simple ‘ lements. By this “package treat- 
ment” of groupings of statical values, proofs and derivations are simplified. The 
force and deformation methods lend themselves to condensation into pairs of 
matrix equations, in which the left-hand sides comprise the stress resultant and deflec- 
tion respec ively. 

For example, derived with the force method, the universally-valid tress resultant 
and deflection matrices, are, as is known, 


S= [J at BY)(BYT EB) -1 BOTFTBOR 


D= SFB 
with 
S,,; — Stress resultant matrix 
D,» — deflection matrix 
Bo) — matrix of stress resultants due to unit redundants 
BY) — matrix of stress resultants due to unit external load 
R,- — load matrix 
F — flexibility matrix 
I  — unit matrix 


where the index ranges are 


a—li.jxnxXp ei— 1... j Sn 
b —1l..i h—ljxnxp 
¢) = 1ij mean Xp Pr —l.j kn Xp 
d —1l..l1xm s —luwk 

e —l..lxm u—l.uuk 

fa TAK vy —L..lxm 
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and i — number of redundants, 
— number of segments, 
— number of loading cases, 
— number of points of application of external loading, 
number of components of external loading, 
— number of contact points of one segment 


Sy 2.9 Bos 
| 


— number of components of resultant stress. 


Ways can be shown for further reduction of the four initial matrices BY, B©, F 
and R into still simpler forms. In connection with the resolution of a structure, 
of continuous rigidity and load, into discrete segments, suitable expressions for 
segment rigidities and nodal loads can be given. Matrices can be partioned for treat- 
ment on small computers. 


| 3. Preparation of flow diagram 
In general, the flow diagram has the following form — 
| 


Input of 

} parameters & data 

| uy 
| 


Setting of indices 


Calc. of address 
| Change of index | 
| 


Logical decisions | 


By the example of a matrix multiplication flow diagram, consisting of a triple loop, 
the concept of programme loop and the role of the index register are explained. The 
matrix multiplication flow diagram shows why matrix operations are especially 
suitable for programming. 

The disadvantages of the matrix method lie in the large memory space required 
and the excessive amount of so-called zero operations (with zero values in the 


Main part of com- 
putation Results 
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memory locations involved) due to the large number of zeros in some of the | 


matrices (esp. the #'-matrix). 


Stability problems can be formulated as matrix eigenvalues. 
Coding, testing and the production are the subsequent stages of the operation. 


Iv. ALTERNATIVES TO MATRIX FORMULATION 


Frequently occurring problems (e.g. continuous beams and frames) can be repre- / 


sented in non-matrix form, in view of the disadvantages referred to above. 


Finally, dimensioning problems involving secondary considerations (¢.g. mini- _ 
mum dead load) can be tackled by means of the decision-making capacity of the 


computer. 


APPENDIX 


As an example of the general form of a flow-diagram mentioned above, a flow 


diagram for a matrix multiplication is described: 


Given: The(m x n) matrix A, 
and the (n x p) matrix B,,; 

Matrix A is stored beginning with memory location R, with successive column 
elements in successive memory locations (411, @12.-. Gm1,@12---)- 
Matrix B is stored beginning with memory location S, according to the same rule 
B54 (Beis B31 > Uy3s-D 49 5-:) 
Required: Matrix C;; = A x B, stored beginning with memory location T, according 
to the above rule(¢p.g: Cin ni, Cis.) 

In mathematical formulation, the ‘‘address functions’’ are expressed as follows: 


>ay& = R+m(k—-1)+(i-1) 
>< = S+n(j —1)+(k-1) 
<= T + mG-1) + G-1) 


where the sign > < means: address of. 


The index register of a computer is capable of adding its contents to the address- 
part of a quantity, and by that to change its address. 


Notation of blocks is as follows (see flow diagram): 


Block I Input of parameters and data 


Block IVa Initialising 
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Blocks II, III, V 
Blocks IV, VI 
Blocks VIII, X, XII 
Blocks IX, XI, XIII 


J. LOEWINGER 


Setting of indices 
Main part of computation 
Logical decisions 


Change of index 
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I Input of matrix elements a;,, bx; 
according to ‘‘address function” 


Input of R, S,T, m,n, p 


Initialise Cij 0 


Form: ¢13<—@y, X byy + ey 


k<ek+1 
Add mto IR2 
Add 1 to IR3 


i<i+1 
Add 1 to JRI = 


3 no 
ps x 


Put m(j—1) + in index register 1 


Put + (i—k) in index register 2 
nj — 1) + in index register 3 


Take a, from memory location (R + contents of IR 2) 
Take b;.; from memory locatioh ($+ contents of JR 3) 


Store c;; at memory location (T + contents of JR 1) 
and keep in accumulator for next cycle. 
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STRESSES AND DEFORMATION 
IN INTERNALLY LOADED CONTAINERS* 


H. ILBERG AND Y. YARNITSKY 
Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


A container internally loaded with a uniform hydrostatic pressure was considered. The 
| theoretically calculated displacements and hoop stresses in the latter were compared 
) with those obtained by means of mechanical measurements, when special attention was 
| paid to the welded seam. 


Because of the wide usage of containers consisting of a cylinder and two hemis- 
_ pherical domes their state of deformation under internal pressure loading has been 


| the subject of intensive experimental as well as theoretical study. This paper treats 


such containers made of SAE steel (Cr, Mo) (Figure 1). 


SE LIZ" 


S 


We e 
POT 6C4S 2 1) O48 34407590" 


1 Figure 1 


_ After flexural rigidity and equilibrium of a spherical shell element is taken into 
) account one obtains the following differential equations: 


dVaatdV. 2 
fo Apia o — (cot*d — pV = SEOR 
a’o dé 2 R 
Sen aa NG ON cen aes aa V 


| * Based on thesis on the same subject by Mr. Y. Yarnitsky, Technion, Haifa, June 1961. 
| Received May 15, 1961. 
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where — 
V — shear force R — radius of shell 
6 — rotation angle of E — Young’s modulus 
shell 
pt — _ Poisson’s ratio D — shell flexural rigidity 
6 — thickness of shell 


Approximate solutions of these can be obtained by using either Geckeler’s! or 
Hetenyi’s! approximation methods. In the first method the functions 8 and V an 
their first derivatives are assumed to be small compared with the second derivatives 
and are accordingly neglected. In the second method only the functions @ and VY 
are neglected. These methods yielded almost identical results when the radial 
displacements were calculated. Neglecting no term in the differential equations: 
a more accurate solution was arrived at by means of numerical relaxation. Th 
differential equations assume the algebraic form: 


Ben re eta — (tan’a — p)Vo = Uo 


Uy + Uz — 2Uo = tana (uy — 


Ur) 
ip A 2’ — (tan’a + pw) Up = CV 


where — 
u = OERO 
R?2 
Cl =s.— D OE 


| = distance between two measured points 


After going through repeated iterations the numerical solution of the differential 
equations was arrived at for the following two cases when the rim of the hemis-: 
pherical dome is loaded by: | 


ins hV ie a0 My, +0 
ble Vo =e 0 Mo = 0 
where: Vo — shearforceattherim M,— momentattherim 


Geckeler’s solutions for the above-mentioned cases were used as an initial approxi- 
mation; this choice proved to be quite successful. Indeed the deviation of Geckeler’s 
solution from that obtained by relaxation indicated that the former was fairly ace. 
curate (Figures 2, 3). 
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2028 -73  -209 +2257 +3246 «257 +1668 ‘76 RELAXATION 
720) -SBh 03.85 3676437 6374 2655 116 4781 #46 216 +186 GECKELER 


pete || | 


7 


hil 


Figure 2a 


1347 0669 0.223 -0.026 -0.132 -0.151 -01275 


Figure 2b 
U and V curves when M ¥ 0, V = 0 according to Geckeler’s and relaxation methods. 
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RELAXATION 
1239 70.062 -11%8 —-1343 “4250 0.858 0524-02 -02564 -0.0832 +0.01266 +20536 0051 GECKELER 


401 Shh 2453 1496 76.22 2905 +415 -516 -5505 RELAXATION 
401 352-2545 1562 624 2675 -2243 1478 -1725 -1465 -1036 Th GECKELER 


belee  L[ie bal Da 


rere wa 
| ea fe 
Rls bitiaeeped 


294 336 378 42 46.2 


42 84 126 168 21 2 


Figure 3 
U and V curves when M = 0, V + 0 out according to Geckeler’s and relaxation methods. 
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Figure 4 


The behaviour of the cylindrical part when loaded by (7) or (ii) was accurately 
calculated assuming that the rims are far apart?. The actual values of Mo and Vo 
were then arrived at by requiring continuity of displacements as well as angle of 
rotation at the cylinder-hemisphere joint. 


Radial displacements were measured over the surface of the container, in particular 
in the vicinity of the weld (which was especially ground in this case) by means of 
an “electronic gauge” (limit of accuracy 1/1000 m/m) with the aid of a device pro- 
viding two planes of reference (Figure 4). 


Measurements were carried out under an internal pressure of 500 atm. The hoop 
stresses calculated from the radial displacement measurements were obtained by 
superposition of the membrane stress and the flexural stress near the weld’. 


Results are given in Figure 5 showing that measured radial displacements can 
serve as a fairly close indication of the hoop stress in the container, in particular 
in the zone of the welded seam. 


For small containers, the method based on measuring the radial displacements 
is preferable to the strain gauge technique for the following reasons: 


|. The order of magnitude of the radial displacements is bigger than that of 
the tangential. Consequently the former can be more easily measured using devices 


which are simpler and more robust than strain gauges. 


54 
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02 
O2m 


Weg 


DY, 


H. ILBERG AND Y. YARNITSKY Bull. Res. Counc. of Israel 


Figure 5 


Stresses and radial displacements in the zone of the welded seam 


— theoretical meridian stresses We, — theoretical radial displacements according 
to Hetenyi 
— theoretical hoop stresses w, — theoretical radial displacements 
— hoop stresses as obtained from the w,, — measured radial displacements 
measurements of w,,, 
theoretical radial displacements according to Geckeler 


Readings of the radial displacement refer to the state of deformation at a point 


and not to the strain integrated along a certain line or area, as would be measured 


by 


Il 
2s 


3: 


a strain gauge. 


REFERENCES 


TIMOSHENKO, S., 1959, Theory of Plates and Shells, McGraw Hill, New York. 
KANTROWITCH, S. B., 1955. Die Festigkeit der Apparate und Maschinen fiir die Chemische Indus- 


trie, Veb Verlag Technik, Berlin 


TIMOSHENKO, S., 1951, Strength of Materials, part If. Van Nostrand, New York, p. 166. 


BEAMS OF VARYING CROSS-SECTION IN BENDING 


| M. MITELMAN 
Division of Mechanics, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


The method described in this paper for calculation of stresses in non-prismatic beams, 
is based on the assumption that the normal stress function o,. is of the third degree and 
depends on the change of the cross-section of the beam. Whence the values of the stresses 
Ty and oa, are derived, as well as the differential equation of the elastic line. 


| 
| 
_ Consider a non-prismatic beam, symmetrical about the x-axis, of rectangular 
_ cross-section (Figure 1), subjected to forces perpendicular to it axis. 


oo 
: 


Figure 1 


0,,-STRESSES 


' The few instances for which solutions derived from the theory of elasticity exist 

(Figure 8 and Figure 9) show that the assumption of planeness of deformed cross- 

sections is not valid here, and the function o, is accordingly assumed to be non- 
linear (Figure 2). 


oe (a, ~ a37;) yo, (1) 
where _M _ 3M. 
ae Te. a Dh? 


‘ dc 
w—coefficient smaller than one, and dependent on Epo 
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Figure 2 


The conditions for determining the coefficients a, and a3 in eq. (1) are: 
a) for y=c we have o, = wo, 


b) by equilibrium fouaa =D [ow dy =4M / 
lA 0 | 


The result obtained is 
0, =4[(5 — 3W)n — S(1 — Wyn? Jo, (2) | 


4 


where 7 = 2. For y = 1, we have a linear distribution of o,.. The maximum of the 


o, stresses for OZ Ws is | 


2 _G-3W)’or | - 5 — 3y 


Omax = 45 (i ae W) 15 (1—W) cx ) Cc 


Ty x-STRESSES 


Considering the equilibrium conditions for a part of the element (Figure 3), we have 


——- 


Te 


Figure 3 
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dN, 
hax (3) 


—Ty = TT. = 


N, = {jo,dA= ss Coe Cah ere) = M 
16 é c Cc 


BA 


whence, by eq. (3) with c a function of x and S = = 
| x 


t= [6 W)~ 25 — 30? + 5 — Wnt | 


— 3 |G - W - 96 - ayn? + 250 - ort | Eo (4) 


a : : pac 
For prismatic beams, with — = 0 and w = 1, the well-known formula for shear- 
ing stresses in rectangular cross-section is obtained from (4) 


3 S 
T tx = 51 — 1’) A 


Verification of the solution according to (4) shows, that [xd A = 1S holds. 
0 


For S = 0, i.e. under pure bending, t,, is not zero. The 7,, stresses at particular 
points of cross-section are: 


3 See a d 
a) y=0 tips =2(5-W)> -s6-WEo 


The contribution of bending to the shearing stress is opposite to that of the shearing 
force. For the loading given in Figure9 1,, =0. 


dc 
b) y=c tye = WO - 


| 
} 
| 
| 
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0,-STRESSES 


Considering the equilibrium condition for a part of the element (Figure 4), we have 


6 ybdx ¥ 


C c+dc 


(5) 


S.= |n.dd=b|td Ltt a(5 Y 4 (5 Layo 30 wy |S. 
vee pra eer yey Oe oe Tele ( at Deis ( rs WaT ( 0 
Hd 


A 


+ Z[G6- W) J 923 we |sem 


whence by (5), with c a function of x and =-q 


oy = ~ 78 — 3S — Vn +15 ~ 30 — 301 — Wyn |Z + 
+ 3 |G - wn -26 - syn + 5 - wy ES 


~ 4 6-Wn—46- ayy + 150-wyns | (Z2)'o, + 


+ §|6- W)n — (5 — 3y)n3 + 5(1 — Wyn l= (6) 
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For prismatic beams we obtain, by eq. (6), the well-known formula for the stresses 
o, in rectangular cross-section 


oy= —4(2— 3-0) 


The a, stresses at particular points are, eq. (6): 


1q dc \2 
a = C28 imeem Se b) y =i-c.. <¢,= (=) 
do! wD G ). y =i + y= + (=) vo, 
2 
q (ce 
c =-C¢ oy= —=-—(—}) a 
ea meh Van) 
PLANE STRESSES FOR y = € 
_ The numerical values of the three stresses (Figure 5) form 
= ($2) v 6- 0, ton’ 
p a 
Txy 6 
aes max 
ad de 6 max 
} 5 3 
| f Tyx ® ge VOr* Ox tana 
| Figure 5 Figure 6 


a geometric series. The principal stresses attained (Figure 6) are: 


F max = (1 + tan2a)Wo,; Omin = 0 


ELASTIC DEFORMATIOM OF THE NEUTRAL AXIS 


By Hooke’s law, we have (Figure 7) 
dp = — —~dx 


Putting y = 0 in eq. (2) after differentiating, we have 


d. 
ab= (5 ~ 30) 
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dp 
&@ 
dv 
dx 
Figure 7 
The differential equation of the elastic line 
d*v 1 M dx 
roles mp deemed cia (7) 


THE NUMBER wW 


An idea of the value of y can be obtained from the solution of the following well- 
known problem in elasticity: 


Ss 
oO 
a 


Figure 8 Figure 9 


a) Wedged cantilever subjected to a moment M (Figure 8) 


Ge 4tan2a tan2« 
~ -3(tan 2a — 2a)(1 + tan2«)? 
b) Wedged cantilever subjected to a load P (Figure 9) 


ve 4tan3a 
~ 3(2a—sin2x) (1 + tan2a)? 
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Results are given in the following table: 


0 


1.000 | 0.9115 0.7889 | 0.7006 | 0.3333 | 0.0800 | 0 


1.0000 | 0.9442 


0.8390 | 0.5833 | 0.3017 | 0 


In order to check the validity of the suggested method the following stresses were 
calculated for the above cases with tana = 4: 
a)o,forn=4 b)t,,forn=0 c)o,forn =. 


Results are compared in the following table: 


Loading according to Figure 8 Loading according to Figure 9 


Exact yy = 0.7 Deviat. Exact Ww =0.839 | Deviat. 


= + 3.8% 3.478 - 3.4235 
2 


4.086 bi bl 


L613 —1.8% 0 0 


Oe i lee +5.0% ra Wes 0.214” 


CONCLUSIONS 


The practical use of the suggested method depends on the correct choice of 
the coefficient w as a function of tan «= aa Taking wy = 1, as is often recom- 


mended, leads to serious deviations in the calculated results. 
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ELASTIC ANALYSIS OF COMPOSITE WALLS - A GENERAL THEORY 


S. ROSENHAUPT 


Structural Division, Building Research Station, Technion-Israel Institute 
of Technology, Haifa 


ABSTRACT 


A general solution for the composite action of masonry walls and the beams or ties stiffen- 
ing their outer edges is given and theoretical solutions for particular cases are derived from 
it. The influence of the position and rigidity ratio of wall and ties is studied. 


INTRODUCTION 


Tests carried out in recent years have shown that steel or reinforced concrete beams 
and stanchions forming the structural framework of a wall do not behave in accor- 
dance with the accepted statical assumptions, and that the filling stiffens the struc- 
ture, reduces deformation under load and increases strength beyond what can be 
expected in elastic or plastic design of the framework. 

This stiffening is of economic and technical importance. From the economic 
viewpoint, it permits reduced dimensioning and saving in the materials required 
for the structural skeleton. From the technical viewpoint, it affects the structural 
behaviour. The beams and stanchions act mainly in tension and compression instead 
of bending, so that the load distribution on the various members as well as the 
distribution of their internal stresses may deviate considerably from the accepted 
assumptions. A study of the composite action is thus of importance for economy 
as well as for closer approximation of the actual behaviour of the structure. 

The study of composite walls began in the last decade and was concerned mainly 
with the practical solution of particular cases, mostly by experimental methods. 
The problems studied were: supporting beams for load-carrying walls,!-4 walls 
subjected to foundation movement 7-12 and to horizontal raking forces within their 
plane, 13-15 and stresses due to temperature and to shrinkage of masonry walls 
restrained against contraction at foundation level 17, Similar problems of composite 
action are also encountered in other fields, but these mainly involve components of 
the same material; for example, in shipbuilding the cooperation of the deck and the 
superstructure 19 and in aeronautical engineering the interaction of stressed skins 
and their reinforcing ribs18. 

The present paper summarises part of the work done on the subject at the Building 
Research Station of the Technion, sponsored by the Housing Department of the 
Israel Ministry of Labour. Its objects were: 
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(a) Development of a general theory of the composite action of walls and the 
beams, ties, and stanchions forming their structural framework. 

(b) Theoretical analysis of a number of problems constituting particular cases 
of the general theory. 

(c) Development of approximate practical methods of calculation based on 
the results obtained by the accurate theoretical solution of particular cases. 

(d) Comparison of theore*ical and experimental results. 

Points (b)-(d) are reported ‘n separate papers, now in preparation. 


STATEMENT OF PROBLEM 


A masonry wall enclosed between reinforced concrete beams, ties or stanchions 
constitutes a load-carrying panel, stiffened at the edges with linear elements and 
subjected to axial forces and bending moments. Cooperation between panel and 
edge members is achieved by shear and normal forces acting at the boundary. 

The fundamental equations of tall beams and of linear elements are known. 
In the case of the composite structure, if boundary conditions between its components 
are written down, and a stress distribution found so as to satisfy both the fundamen- 
‘al equation and boundary conditions—a general method will be provided for the 
solution of a series of problems hitherto treated separately but which are in fact 
particular cases of the general solution (such as load distribution in foundation 
beams, differential settlement etc). The loaded structure will be represented by the 
fundamental equations and boundary conditions, and the particular case by the load 
or external deformations imposed on it. 


METHOD OF SOLUTION 
Fundamental equations 


Panel. The Airy function for the panel is defined as the function whose second partial 
derivatives are the internal stresses in the panel: 


ee a aoe 020 
Ox = “Gy? PIs Oe x =~ Oxdy 
the fundamental equation for the panel being: 
0+® 04+ 04 
—— —s +z = 0 
Ox4 +2 0x20 y2 - Oy* 
Beam. The fundamental equation of a beam in bending is 
02v 
M, ee E,I, ae 


where E, = modulus of elasticity 
I, = moment of inertia 
M, = bending moment 
= displacement normal to beam axis. 


Ss 
i 
| 
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Tie. The equation of a tie subjected to an axial force is: 
Ou, 
P; = Ag 40x = A,E, ins 
where P, = axial force 
A, = beam cross-sectional area 


u, = displacement in direction of tie axis. 


In the following analysis Poisson’s ratio will be neglected (v =0). Actually, 
its inclusion would not complicate the process, but it was found that such a degree of 
accuracy would not be required, considering that of the basic data. 


Boundary conditions 


The object of the following derivation is formulation of the boundary conditions 
between wall and beam as an expression containing the Airy function (®) of the 
wall only, so as to deal with only one unknown at every interior or boundary point 
of the panel. For this purpose, the conditions of equilibrium and displacement at 
the boundary will be written down; then, by the laws of elasticity relating forces 
and displacements, elimination of the former would give the boundary conditions 
solely in terms of displacements. The latter, in turn, will be expressed as functions 
of ®, thus giving the boundary conditions in the required form. 


YY 
we 
9) 


Figure 1 
Forces and moments acting on the beam at junction of beam and panel 


Notation (Figure 1) 
P, — axial force in edge beam 
M, — bending momentin edge beam 
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N, — tension per unitlength between paneland edge beam 
Ny, — shearing force per unit length between panel and edge beam 
p  — external distributed force applied to beam directly 
S — shearing forcein edge beam 
u — displacement parallel to direction of edge element (x) 
v — displacement normal to direction of edge element (y) 
subscript s — referring to edge of panel 
subscript b — referring to centre line of beam 
t — thickness of panel 
A, — Cross-sectional area of beam 
d, — depth of edge-beam 
£,,€, — fibre strains in x-,y-directions respectively 
GO, 0, — Stresses in x-, y-directions respectively 
Txy — Shearing stress 
Yxy — Shearing strain 


Equilibrium conditions 


The equilibrium conditions of an element of the edge beam (Figure 1) are: 


inx -direction: dP — N,,dx =0 (1) 
in y-direction: pdx + dS — N,dx = 0 (2) 
d, 
moments about z-axis: dM, +N,,dx ra — Sdx = 0 (3) 


Conditions (1), (2) and (3) hold for elastic as well as non-elastic materials and for 
both small and large displacements. 


Displacements 
The conditions of common displacement of the panel and edge beam at the junction 
(Figure 2) are given by: . sR, ; 
pir a ae oe (4) 
i =, (5) 


Conditions (4) and (5) hold for smalidisplacements only. 


Figure 2 
Compatibility of displacements at junction of beam and panel 
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Elasticity conditions 
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Cup ou 
P=Ay5,0, = A,E, ao since o, = Ee, and ¢, = aa 
020 
My sire 


where ,0, = tensile stress at centre of beam. 


Nyx = —ttay = eck ee ge ee Ries so 
eee ie. OV. KOxeer Des Set DANO Vs: SOX 
Ov, 
N, = toy = ae 


Boundary conditions expressed as functions of u, and v 


From eq.(1) we obtain 


dP, 
~ NY = ae 
by eqs. (6), (4) and (5) 
OV ery | Cun. CAE, Ou, _ 02 ov, d 
de da Be ets a ee, AO Bee | ce 3) 


and by eqs. (8) and (1’) 


Omitting the subscript ‘‘s’’ and rearranging: 


LES | Cus —Y CO\ 82 02 
2 a as as ) = Asks ase 


. ou OD _ 2A,E, { 0?u_ dy O3v \_ 
_ Ox? 2 Ox? ) 


‘éy Ox tE, 
b (3 
y eq. (3) 5 LaMegmwapileni anny 
~ ax pe aay 


Differentiating eq. (3), and by eq. (2), (9), (6) and (7) 


dS _@M, , dy dP, 


=—- =Ny-p= Ov, 


= tE 


dy dP, 
DQ ape 


Replacing “‘b”’ with ‘‘s’’ and referring to eq. (4) and (5): 


04v, dy 


dui ici ails ie ak ca as aa 


A,E, 


(6) 
(1) 


(10) 


O3u, 


0x3 
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Again omitting the subscript ‘‘s’’ 


Ov at DSU Ady 03 Ov dy 
Pan p= EJ ea ae 2D A,E, Axe (« = ax. ) 
And finally, dividing by tE, and rearranging: 
a p _a&A pa av (dy | 2 Sy (11) 
OVGgathi tts 2.1) Oxo eax \., ded; 


Eqs. (10) and (11) are the boundary conditions expressed as functions of the dis- 
placements of the panel. 


The next step, as explained earlier, will be to express the displacements in terms 
of the Airy function ®. 
Displacements u, v and their derivatives as functions of ® 


Eqs. (10) and (11) contain the following derivatives: 


Ou dv, 02u_ Ou. dv. O30, O40 
dy Ox’ Ox?’ dx3’ dy’ ax3’ ax4 


which will now be expressed in terms of ® by successive introduction of the functions 
relating ¢, o, and O. 


Ou Ov Tah ee 2 020 

— —— = -—— x = — = — — 12 

oy ea dx 2G E E 0xdy “ 
on a ou) 0 0 6G. Aide, 000 1 03® 


-2 45. = 2 (13) 
xe OX Ox.” OG EE cx E 0x dy? E daxdy? 


us ay 04+@ (14) 
ax3—E axdy? 
C0 tape. 4, ol ee (15) 
ay. 7 OO Le CEES Ox? 
ov ou 2 ou 


a)? ay EE” ay 
C202 BOT Ly 02u ZOU (Oty aed OT, iN o 


ax2—«=E Ax dxdy Ox oy E. Ox bs Oy, 


nN (Rete Os 
E Ox oy 
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2 (16) 


Ae! > O20 - 0 Sn a\ ew ke > eo 0 = 
her dx2 Oxdy}) E Ox3dy dxdy3 
04 04+® 


ee ( reine snaye | 


4 5 5@ | 
A ee Ccree ) (17) 


Ox4 E \~ dx4dy ° dx2dy3 


Boundary conditions as function of ® 


Substituting (12) to (17) in (10) and (11), the former becomes: 


2) TOMD TP 2ALE, pts we Oy diet Os ED 
E.0x0y ~ tE, “VE, “dxdy2d, 228, Ox3dy = dxdy3 


and after simplifying and rearranging, we obtain 


(10 ) 


4 Ox3dy <dxdy3 


LE ee Oy  e2@ d, 04+® 04® 
A,E, Oxdy  dxdy2 2 


Eq. (11) similarly becomes 


EROxt? “4ES 0 iB. 


eure i 


A 02@ = Pp d, A,E, tes ona) aL 0>® + EXO) ) d, cal 
E, Ox2dy? E, ( Ox40y = dx2dy3 ( 


and after simplifying and rearranging, we obtain 
FES 24) 02D ip | et ds a7 a5@ a5@ al 
ArEs dy E : | ~ Ox26y? (3 fi ids | Pees i sna) a 


Numerical solution 


The method adopted for numerical solution was that of finite differences. This 
permits expression of the boundary conditions in all directions, with the degree of 
accuracy of the same order as that of the fundamental data of the problem. 
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Boundary conditions expressed in terms of finite differences 
The derivatives appearing in the boundary condition equations are: 
20 70 80 40 40 4b OO ASO 
0x2” Oxdy’ Oxdy?’ Axdy3’ Ax3dy’ Ax2dy2’ Ax2dy3’? ax4dy 
Using the coordinate grid shown in Figure 3, with its origin at 0 and h and k as 
spacing units, the above derivatives are determined as follows: 
Co @, — 20-0, 
Bry maa ue) 


020 oO, + D, oan 0; — @, 
3 ee eee 19 
axdy 4hk (1) 


o3® ts) 02@ ®, —— 20, + ®, _— (®, aay 20, + ®,) 
CN Ge = 20) 
Oxdy2 dx \ dy2 2hk 
_ 0, + ® — ® — O, — AD, — O,) 
- 2hk2 
@xdy3 ox \ dy3 ] ~ 2h-2k3 
_ Dy, + Dyg — D5 — Dig — APs + O, — De — De 
2 4hk3 
Similarly: . 
AD D1, + O17 — Di5 — Dao — 20, + ®, — O, — G3) (22) 
Ox3dy iz 4h3k 
04D 02 (02D\ OD, + Og — 20; — 2, + O, — 2H) + 0, +0, — 20, “i 
Ox2dy2  Ox2\ dy? F h2k2 
_ 4@y — 2, + ®, + D3 + @,) + (M5 + Oe + O, + Dp) (23) 


h2k2 
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5b a2 (03D \ _ ee 
@x2dy3 ax? \ dy3 J 


O14 — 205 + 2Dg—Dy 9 — Ayo — 2D, + 2H, — Oy,) + Oys— 205 +20,— Dig _ 
h2 + 2k3 


Di, + D5 — Dig — Do — 2@; + O, — O, — O, + Di, — O,,) — 4, — OG.) 


2h2k3 
a _ 0 (a4) _ 
axtdy dy \dx4 ] (25) 
_ O13 — 405 + 6@, — 406 + Dig — (Mr — 405 + 6, — 40, + O17) = 
= 2kh4 


®13 + Dig — Diz — Pro — 4s + Oe — O, — Dy) + OM, — O,) 
2hk4 


Figure 3 


General coordinates for finite difference approximations 


Finally, substituting (18) to (25) in eqs. (10’) and (11’), the boundary conditions 
are obtained in terms of finite differences 


ktE, 
24, Es (®,; ci 0, 7% OD, ma ®,) apa (®, a @, Sail D, on @, = 20, = 20,) + (10) 


dyk 
+ are + D7 — 1, — ®,) — Ab; + ©, — O — O,)] 


d 
us az [®,, + Di, — D5 — Dy — 210, + ©, — ©, - ®,)| 
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KE, 2 h2 
AsEy dy (®: - 2% +0, - pF) = a 


= 4o — 2, + O, + ©; + D4) + (OD; + O, + O, + Oy) + 


dae ny te 
AG = id; iis [®1s + ®y¢ — Py — Oyo — 4(@, + ©, — ©, — @,) + 


if 


+ 6(®, — ®,)| i [Ors — 91, + ©; — O,, + 4(@, — ®,) — : 


— ®, + O, — O, — O; + D1) — ®,.)]} 


SIMPLIFIED FORM OF SOLUTION 


Kgs. (10) and (11), in all their various forms, are the complete mathematical expres- 
sion of the composite action. These equations have been derived and written down 
in full, so that, if required, the effect of the intuitive disregard of some of the terms 
could be easily ascertained. 


In their complete form, the equations are obviously too complicated for practical 
use. A justifiable simplification would seem to be to neglect the rigidity in bending 
of the beam compared with that of the wall, which 1s 8 to 10 times higher than the 
beam, so that its moment of inertia is 500-1000 times larger. In these circumstances, 
it cannot be expected that the resistance of the beam to bending will have a consi- 
derable effect on the internal stresses of the wall. Disregard of this resistance is ex- 
pressed in the general equations by equating the moment of inertia and depth of the 
beam to zero. 


This is done as follows: 
Assuming d, = 0 and J, = 0 in the boundary condition (10’) and (11’), we obtain: 


tE, 020 3® 
aercaSs = 10’” 
A,E, Oxdy  Oxdy? oo) 
020 Pp 

They ct 0 1122 
Ox? t ( ) 


The above equations can also be obtained directly (Figure 4), as follows: 
Equilibrium conditions: 


pdx — N,dx = 0 (27) 
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but: 


020 
tS =f 28 
Nyx ITxy t Oxy ( ) 
and: 
E, 
P, = Ay 5x a Ay & 3% (29) 


Differentiating equation (29) with respect to x: 


OP, _ AsE, 0,0, AE, 0 (@\_ A,E, 030 (30) 
Ox iy Ee. Ox Ee Ox Ny?) Ee dxdy? 
Since 
62@ 
sox Oy? (31) 


Substituting eqs. (28) and (30) in eq. (26): 


E,A, 20 _ j 020 
E, Oxdy2  @xdy 


OF; 
, 80 33 


A,E, Oxdy  axdy2 


The above is identical with condition (10’’’) obtained earlier as a particular case of 
the general equation. On the other hand: 
020 


N, = to, =t Ox2 (32) 


Substituting eq. (32) in eq. (27) we obtain 


= 020 
Ox? 


P =) 


This in turn, is identical with condition (11”) similarly derived earlier as a particular 
case of the general equation. 


Eqs. (10’”) and (11’’) can be expressed in terms of finite differences by using 
formulas (18), (19), (20), as follows: 


ktE, 
2A,E, (8 + ®, — ® — ®) = 0, + ® — ® — O, — 20, + 20, (10°Y) 


a, cia 2D5 = 0, Pp 
—ecre (11%) 


| 
| 
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Parameters 


A group of coefficients appearing in all fundamental equations (10) and (11) in 
their various forms is: 


tE, 
B na Arr. (33) 


This parameter relates the rigidities of the panel and edge beam against normal 
forces of tension or compression. As the rigidity of the beam increases, 6 decreases. 


Particular cases: 


1) Edge beam of infinite rigidity: A,E, > 0 £0 
2) No edge beam: AE, = 0 Bo 


A non-dimensional parameter «, given by: 


a 


1 vEy 


o> ag Sadh. 


(d = height of wall panel), (33’) 


appears in all numerical calculations as wellasin the proposed approximate methods. 


| PARTICULAR CASES. 

- On the basis of the general theoretical solution derived in the previous sections 
a number of particular cases was studied. Two of them are summarised in the 
following. 


a) Load-carrying wall on point-supported foundation beam 
_ This distribution of the load transmitted by a wall to its supporting beam is material- 
_ ly altered by the “‘arch effect”’ taking place in the wall. Among the methods proposed 
_ for the practical solution of the problem may be cited the triangular load diagram 
- method’, an empirical method of moment coefficients* and a theoretical method 
- based on the radial diffusion of the column load in the wall (considered as a semi- 
infinite plane)‘. 
In the present study wall and beam were regarded as components of a composite 
structure connected by the boundary conditions as derived in the preceding para- 
graphs, the case of the load-carrying wall on supporting beams constituting a parti- 
cular case of the general solution. 


The cases studied were: 


Height-to-length ratio of wall 2:3 (corresponding to that of laboratory models) 
with modular ratios E,/E, of wall and beam as follows: 


74 S. ROSENHAUPT Bull. Res. Counc. of Israel 


1) & (Ytong block wall on reinforced concrete beam); 


>) 535 “ (burnt brick wall on reinforced concrete beam) ; 


3) 1 theoretical ratio, for a more complete picture of the influence of the 


modular ratio). 


The values of « (formula (33)) corresponding to the above cases are respectively: | 


(1) 3.580; (2) 0.628; (3) 0.358. 


The case of an homogeneous wall (a = 0) was solved for comparison. The dis- | 
tribution of o,, stresses along the central section of wall for different beam-wall | 
rigidity ratios (a) is given in Figure 4. It is seen that as beam rigidity increases, the 


tensile stresses at the bottom of the masonry decrease, as all tension tends to con- 
centrate in the foundation beam. 


Figure 4 


Forces acting on a tie at junction of tie and panel 


Point-supported wall with lateral ties. 


In practice, walls strengthened along their lateral edges are used in prismatic water _ 
towers and silos. In such structures, the wall is usually regarded as a diaphragm | 
girder. In its stress analysis the contribution of the perpendicular walls is normally | 
neglected. It is to be supposed, however, that the o, stresses are concentrated at 


the rigid corners of the wall, and that this leads to redistribution of o, (which deter- 


mines the main reinforcement) and t,, (which determines the shear reinforcement) | 


in the wall. 
Three numerical examples were solved: 


a) No vertical tie («=0). b) Vertical tie of section area equal to 1/3 of wall’s hori- 
zontal section (~=1/3). c) Vertical tie of infinite rigidity (« =o). 


| 
| 
| 
| 
| 
| 
| 
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The distribution of o, stresses at midspan is given in Figure 6. It is seen that the 
existence of a lateral tie results in considerably lower tensile stresses at the bottom 
of the wall. Instead of concentrating at the bottom, the tensile stresses are distribu- 
ted nearly throughout the entire height of the wall. This explains, in a way, the 
cracks appearing in silo walls in the central zone where membrane stresses are added 
to the tension due to horizontal bending stresses. These membrane stresses are not 
considered in statical computations, which take into account the action of the com- 
mon diaphragm girder (and not a composite girder as existing in reality). 


-28 -24 ~=20 


-—— COMPRESSION TENSION ——~— 
Notation: b 
—__—_—— Gh (0) fn 
gg a= 0.358 _ A,E, Es 
eer nnn ene a = 0.628 oe tdE ? i 

G7==)3580 © a 
Figure 5 3 


Wall on supporting beam. o,, stresses in midspan section expressed 

as coefficients of p/t. Influence of variation of beam-wall rigidity ratio a. ; 

fay |e 
b 


Figure 5a 


It may be observed that the same problem of lateral ties has been studied by 
Gruning 24 by a method of superposition requiring the solution of nine particular 
cases of loading. The results obtained in both studies were compared and agreement 
is good. 
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+40 -08 -06 -04°-02 0 02,04 O06 O8 0 t2 14 16 18 20 
=—— COMPRESSION TENSION ——= 
Notation: 
a=0 = A,E, 
sesasenanacnaaeanese Q = 1/3 tdE 
L= iS 
Figure 6 


Wall with lateral ties. o,, stresses in midspan section A-A 
expressed as coefficients of p/t. Influence of variation of 


Ne 


se 


So OST Tiga ot 


tie-wall rigidity ratio a. 


a=o | section B-B 


Figure 6a 
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A NUMERICAL METHOD FOR THE DETERMINATION OF 
STRESSES IN WALLS WITH RECTANGULAR OPENINGS 


G. MULLER 


Structural Division, Building Research Station, Technion — Israel Institute of Technology, 
Haifa 


ABSTRACT 


A practical numerical method for the determination of stresses in walls is presented. Thd 
method was developed for the determination of stresses in walls with openings. It ca 
however be adapted to a great many biharmonic problems involving internal boundaries 


1. DEFINITION OF THF PROBLEM 


Stresses in walls are governed by the Airy stress function ¢ which satisfies th 
biharmonic equations 


ot 4 4 
do , oe 
Ox4 Ox2dy2— 0y4 


The stress components are 


6, = 0d , 
Oy? () 
oat | 

C153 a (2 
ad 

pa Oxdy (3 


In order to calculate the values of # at any point on the surface of the wall, i 
is necessary to define the boundary conditions resulting from the external loading. 
The problem of walls with openings is complicated by the fact that there exist two! 
boundaries: an external and an internal one. The boundary conditions at the external 
boundary can be obtained by integration from the known forces distributed along 
this boundary, and the constants of integration may be chosen arbitrarily. But once: 
these constants of integration at the external boundary have been determined, the: 
corresponding constants at the internal boundary cannot be chosen in the samc! 
manner, and before the function @ can be defined, it is necessary to determine the: 
constants of integration at the internal boundary. 
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2. THE BOUNDARY CONDITIONS 


The boundary conditions at the external boundaries depend on the external forces 
jand can be determined by any of the known methods given by Timoshenkol, 
Girkmann?, Rosenhaupt3 and others. In this paper, the boundary conditions at 
the internal boundary only will be discussed. 

The boundary conditions at the internal boundaries are obtained from the condi- 
tion that at the opening, shear and normal stresses are zero. Taking a rectangular 
|Opening, 7, and t vanish at the vertical edges; o, and t vanish at the horizontal 
sedges. # is obtained by double integration of tie stress components, and in order 
{to satisfy the above conditions, ¢ must vary linearly along this boundary. 

The boundary conditions at the internal boundary are therefore: 


o = Ax + By +C (4) 

0 

aed ©) 
| eu? (6 


)A, B and C being constants of integration which are still to be determined. 


ig 


3. OUTLINE OF PROCEDURE 
'The values of ¢ are calculated by the finite difference method, and the constants of 


‘integration A, Band C at the internal boundary are obtained from the conditions of 
statical equilibrium. 

The problem is thus reduced to the solution of two sets of linear algebraic equa- 
tions. On the other hand, the stress function ¢ ,instead of being satisfied at any 
point on the surface, will be found at the nodal points of the finite-difference net 
‘only, and the accuracy of the solution will therefore depend on the fineness of the 
net chosen. 

In Figure 1 a square finite-difference net is given. The nodal points are numbered 
from (0) to (12). The spacing of the nodes is h; and dp .... 12 denote the values 
of # at the corresponding nodal points. 

The Airy stress function, expressed by finite differences, becomes an algebraic 


equation: 


20b0-8(b+$2+b3+ bs) +2b5+ be + $7 +¢s)+ bot biot+ $11+¢12=0 (7) 


By writing condition (7) for every nodal point on the wall and substituting 
relevant boundary values, a set of simultaneous algebraic equations for the calculation 
of ¢ is obtained. 
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42 


Figure 1 
General coordinates for finite difference approximations 


In walls with openings, the boundary values at the internal boundary are not defined.. 
They are given in terms of the constants of integration A, B and C, which are still] 


to be determined. The complete solution of @ cannot therefore be obtained unless} 
these constants are known. It is however possible to obtain partial solutions for’ 
go, O4, Pg and ¢¢ from which partial stress components can be computed. 
do is the value of ¢ due to the external loading P, ¢, is the influence value | 
of ¢ due to the integration constant A, and ¢, and ¢¢ are the influence values of @ | 
due to B and C respectively. 


These partial solutions are obtained by solving the following four cases: 


(a) For¢?o: P=1 A=+O B=0 and C=0 
(b) For¢@,: P=0 A=1 B=0 and C= 0 
(c) For ¢,: -P='0 “A= 09 BEd tandy C0 
(d) For¢?e: P=0.A=0 B=0 and C=1 
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The complete solution of ¢ will be 


} = $o + Ady + Bhg + Che, (8) 


A, B and C are still to be found. 
The stress components are obtained analogously. 


Cx x0 am AOx4 23 Bo, 7 Cox (9) 
Oy => Oyo J. Ady4 ap Boyc + Coyc (10) 
) T=) + Aty + Btz + Cte (11) 


© 50> Ox» OxB and o,¢ etc. are the partial stress components due to ¢o, by, bg and dc 
respectively. 


4. THE CALCULATION OF A, B AND C. 


The constants of integration A, B and C are obtained from the conditions of 


' equilibrium. These conditions are: 
y Mi= 04>. V = 0 andy = 0 


For any vertical cross section, these conditions become: 


YM = Yo,y(Ay) + M. =0 (12) 
! LV = YrAy) + Q, =0 (13) 
| LH = Yo,(Ay) + Q,=0 (14) 


| where M,, Q, and Q, denote the external moment, vertical shear and horizontal 
shear, and (Ay) denotes the vertical spacing of the nodes. These conditions may 
| be written in the following form: 


| Y opoWAy) + AY oxsy(Ay) + BY. oxgy(Ay) + CY oxycy(Ay) + M. = 0 (15) 


> to(Ay) + AY ta(Ay) + BY te(Ay) + CY tAy) + Q, = 0 (16) 
EY, o,o(Ay) + AY, o,4(Ay) + BY) oxx(Ay) + CD, oxc(Ay) + Q. = 0 (17) 


It should be noted that all values required for the formulation of the above condi- 
tions must be calculated by the trapezoidal rule. 

Conditions (15), (16) and (17) are three equations containing A, B and C, from 
which these constants can be determined. Furthermore, as these conditions hold for 
any vertical cross section, it is possible to solve walls with any number of openings. 
There are three constants ofintegration for each opening, and the number of equations 
required will be three times the number of openings. 


| 
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5. FINITE DIFFERENCE APPROXIMATIONS | 
For the formulation of conditions (15), (16) and (17) and for setting up the boundary ff 
conditions, a number of finite difference expressions is required. A list of these) 
expressions is given below. Notations refer to Figure 1. 


dp od: — $3 
est Voy one eee ae, 18 
(5 ), 2h | 
06) _ b2 — ba | 
eh iat i eI 19)» 
( oy ), 2h co | 
0° _ 1 + 63 — 20 \| 
Cae iF a h2 (20)| 
arg 2 + 4 — 2b | 
ae ye 2 (21) 
ao $s — bo + $1 — Ys . 
far iz 4h2 (22) 
6. THE BOUNDARY VALUES IN FINITE DIFFERENCE FORM 
At the internal boundary, ¢ = Ax + By + C. In order to express = and at at | 


the boundary by finite differences, imaginary values of ¢@ inside the opening must | 
be introduced. | 
Let ¢,, be the real value of ¢ at the wall near the vertical edge of the opening, | 
and ¢,, the corresponding imaginary value of ¢ inside the opening. 
dp Pev = Piv dp 


By (18), ae SEARLS But ea A, hence ¢;, = ¢., — 2Ah, and | 


ifh=1, dy, = Vee 2A 


Similarly let @,, be the real value of ¢@ at the wall near the horizontal edge of | 
the opening, and ¢;, the corresponding imaginary value of ¢ inside the opening. 


Op _ ben — Pin Ct | 
By (19), Apia Weale as aye B, and 


bin = en — 2Ah = be, — 2A. | 


A typical scheme of the boundary values at theinternal boundary is given in Figure 2. | 
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Figure 2 
Boundary values at internal boundary of wall with opening 
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7. SINGULAR POINTS | 

The corners of the opening constitute singular points at which the boundary 

conditions cannot be strictly satisfied. A certain error is thus introduced into the 

calculations, which makes it impossible to compute the stresses in the immediate | 
vicinity of the corners. The magnitude of the error depends on the finite difference 
net chosen. By working on a sufficiently fine net, i.e. keeping / sufficiently small, 


any required accuracy can in fact be obtained. 


By the principle of St. Venant, the error will not affect the stresses at some distance 


from the corners. The constants of integration will therefore not be affected if the 


vertical cross section at which they are determined is sufficiently removed from the | 


corners of the opening. 


8. SOLUTION BY MATRIX INVERSION 


Walls with openings require a very much finer finite difference net than walls without 
openings. It is often required to solve sets of 30 or more simultaneous linear equations. 
It is therefore convenient to solve the equations by matrix inversion. 

Let the stiffness matrix of the wall be denoted by [A]. This matrix depends on the 
shape of the wall alone, and is independent of the boundary values. The boundary 
values are denoted by the matrix [k]. The equations are then given in the form 
[A] [6] = [k] and [¢] = [A~!] [kK] where [A~!] is the inverse of the stiffness matix 
[A]. The advantage of this method is that once [A~1] has been calculated, @ for any 
given boundary conditions [k] can easily be obtained by simple matrix multiplication*. 
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* The method is explained in detail in P. B. Morice’s ‘Linear Structural Analysis’. 


STRENGTH AND DURABILITY OF MORTARS MADE WITH 
AN ADMIXTURE OF LOCAL POZZUOLANA 


R. MALINOWSKI 
Building Research Station, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 
Technical properties of Carmel pozzuolana as admixture to concretes and mortars for 
hydraulic structures were studied. Uncalcinated pozzuolana was ground to a fineness of 
4800 cm2/g. An admixture of 15% by weight —calculated on cement basis —in 1:3 mortar 
(cement : sand) produced an increase in flexural and compressive strength and improved 
the impermeability of concrete. The increase in sulphate resistance is not found in the 
pozzuolana-containing mortars from the sources studied (uncalcinated). Conclusions on 
durability testing are given. 


Durability of the concrete is a prime factor in maritime and hydraulic structures. 

The present study was undertaken following the discovery of the high durability 
of ancient hydraulic structures in the vicinity of Caesarea and Tiberias (in which 
pozzuolana and other durability-improving admixtures had been used) — on the 
one hand, and in view of the rapid development of local maritime and hydraulic 
construction on the other, as well as on grounds of utilisation of cheap local raw 
material. 

Mechanical properties, impermeability and durability of mortars made with 
admixture of local (Carmel) pozzuolana were studied. Mortars and concretes made 
with 15% and 30% pozzuolana replacing cement (calc. on cement basis) were 
compared in laboratory tests with those made with cement only (ordinary and 
sulphate-resistant). 

The physical and chemical properties of pozzuolana were studied. The pozzuola- 
na was ground to a fineness of 4800 cm2/g (Blaine test). The mortars were cured 
under identical conditions —1 day in air and 6 days in — water and afterwards in 
artificial sulate solphutions: MgSO, (5%) and Na SO, (10%). 

The aspects studied were — 

a) Influence on compressive and flexural strength of mortars. 
b) Resistance of mortars to sulphate attack (mechanical tests). 
c) Water impermeability of concretes. 

Recommendations on the utilisation of pozzuolana and on the direction of future 
research were made, and a durability test (based on a combination of visual exami- 
nation and compression and flexure tests) was proposed. 


1. INFLUENCE ON STRENGTH 
Under water-curing conditions, an admixture of 15% pozzuolana in a 1:3 (cement: 
sand) mortar produced an increase in flexural and especially in compressive strength. 
An admixture of 15° quartz sand also increased the strength of the mortar. 
Results do not permit separation of the chemical and physical effects, so that the 
chemical activity of the pozzuolana is difficult to evaluate. Continuous production 
control in utilising pozzuolana as a substitute for cement is desirable. 
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2. RESISTANCE TO SULPHATE ATTACK | 
On the basis of the three criteria (visual changes and changes in flexural and compres- 
sive strength) the following conclusions can be drawn: 

A 15% admixture of pozzuolana gave positive results on 150 days’ storage — no | 
reduction in strength and considerably superior durability compared with ordinary 
Portland cement mortar. Similar results were obtained with quartz sand. 

On 240 days’ storage a considerable decrease in durability was observed in mortars 
made both with pozzuolana and with ordinary “Nesher’ cement. Mortars _ 
made with 15°% quartz sand and with special cement gave satisfactory results. An _ 
increase in pozzoulana content above 15% resulted in further decrease in durability. | 


3. WATER IMPERMEABILITY OF CONCRETES 
Improved impermeability was observed in concretes with 250 kg/m> and 300 kg/m3 
cement contents on admixture of 15% pozzuolana. Ordinary 300 kg/m? concrete | 
showed water penetration at 3 atm pressure, while no penetration was observed in | 
the admixed concrete at 5 atm. 


4. CONCLUSIONS ON THE UTILISATION OF POZZUOLANA 
Results show that local pozzuolana can be used, without calcination, only as an. 
impermeability-improving substitute for cement in mortars and concretes, in amounts 
of 15%. This conclusion should be systematically verified under pilot-plant produc- | 
tion conditions, e.g. on blocks and other concrete products as well as on wall conc- 
rete made with rubble aggregate and on impermeable concretes for hydraulic 
structures. | 
In spite of the satisfactory behaviour observed under rigorous conditions of’ 
sulphate attack after 90 and 150 days’ storage, no positive results could be obtained 
with the material tested on further storage, so that untreated pozzuolana from the | 
source in question cannot be recommended as a durability-improving admixture. | 
Further study is in order (especially as to sulphate resistance) on material from. 
deeper excavations, as wellas after calcination. Other sources should also be examined. | 


5. CONCLUSIONS ON DURABILITY TESTING OF MORTARS AND CONCRETES SUBJECTED TO | 

SULPHATE ATTACK | 
Evaluation of the durability properties cannot be based on a single factor (such as | 
compressive strength). Surface changes have little influence on strength in view of 
the small thickness of the affected layer. 

Visual examination as well as study of surface softening (incl. hardness tests) 
provides good control. The same is true of the flexural test, in view of the influence | 
of the affected layer on the section modulus. Expansion tests are unsatisfactory in| 
view of the smallness of the deformations involved compared with the dimensions 
of the specimens ordinarily used. 


A complete indicator is provided by a combination of all three criteria. 
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THE RECOVERY OF CRIMP OF WOOL FIBRES 


J. BANBAJI, E. ALEXANDER AND M. LEWIN 
The Institute for Fibres and Forest Products Research, Ministry of Commerce and 


Industry, Jerusalem 


SUMMARY 


The load extension curves of fibres are widely used to illustrate the different perfor- 
mances of various fibres under mechanical or any other treatment. Generally, any 
modification in the inherent property of a fibre is reflected in the load-extension 
curve. The changes in the shape of this curve, under the influence of various treat- 
ments, are in most cases, an expression of variations in parameters such as Young’s 
Modulus, extension etc. 


In the wool fibre, the load-extension curve is composed of three successive regions. 
The first is the initial uncrimping energy which corresponds to the visual crimp 
of the fibre. This is followed by the linear or Hookean region. Finally, there is the 
plastic region which extends to the breaking point of the fibre. The behaviour 
of the fibres in the last two stages is usually explained as being the result of changes 
in the molecular bonds. 

The importance of the uncrimping region is especially great in the course of the 
processsing of wool. In previous works, no attention was paid, for example, to the 
recovery of the crimp and its relation to other properties of the fibres. Also, in this 
interval of low forces applied on the fibre, no attempt was made to analyze the curve 
and to interpret the nature of the forces participating in the process of uncrimping 
extension. The reason seems to be the lack of the required sensitivity of the instru- 
ments used on previous occasions, which is necessary in order to study the curve 
in the low load interval. On the other hand, the load-extension curves of wool 
at applied forces higher than those applied in the uncrimping region, have been 
the subject of an extremely wide range of experiments, which have provided a lucid 
interpretation of the elastic properties of the fibre. It should be noted that in pre- 
vious studies of the recovery of the stress-strain curves of wool fibres, all the expe- 
riments undertaken related only to the elastic and length recoveries. 


* This work is a part of a research project performed under a grant of the International Wool 
Secretariat, London. Full details of this work will be sent for publication to the “Textile Research 


Journal’, U.S.A. 
Received May 15, 1961. 
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The purpose of this work is to investigate the uncrimping region from the lcad- _ 


extension cuive and to study the recovery and stability of the uncrimping energy 


in successive cycles, after removing the load, as related to the percent of extension — 


in the fibre. 

The uncrimping energy is defined as the energy necessary to straighten the fibre. 
It is measured as the initial part of the load extension curve. 

At extension below the yield point, which occurs at 2% of the uncrimped length 
of the fibre, it was found that the uncrimping energy already decreases in the second 


cycle which is performed immediately after the first. For additional cycles there 


is a tendency for a slight and irregular decrease. A typical decrease of the uncrimping 
energy of one fibre 1s given as follows: 


Cycle Ist 2nd Sth 10th 20th 25th 30th 
Uncr. Energy 100.0 91.6 96.2 88.6 92.0 85.0 84.7 


The initial value of the uncrimping energy, when measured on 10 fibres, was 
reduced in the second cycle by up to 15%. 

Beyond the yield point, and for extensions between 2% and 30% of the uncrimped 
length of the fibre, the uncrimping energy decreases after the first cycle, at a rate 
which depends on the percent of extension applied on the fibre. At an extension 
of 30%, i.e. slightly before the breaking of the fibre, it was found that most of the 
uncrimping energy disappears in the second cycle. 


The above experiments were undertaken on two differently crimped types of | 
Merino wool fibres. The results showed similar tendencies for the decreasing of © 


the uncrimping energies on both samples. The following table shows the r2sults 
obtained on one type of fibres: 


Extension of the Recovered Coefficient 
fibre (%) Uncrimping Energy of Variation 

(%) (nd cycle) Yo 

2 85.9 13.7 

5 61.0 14.4 

10 39.7 122 

20 18.0 DPA 

30 6.9 58.0 


This phenomenon was also examined for various times of recovery between the 
Ist and 2nd cycle. It was shown that after one hour of relaxation of th: fibre, a 
greater part of the uncrimping energy is recovered than in the case where the second 
cycle was performed immediately. It was found that at an extension of SO, of 
the length of the fibre, the uncrimping energy, recovers up to 70% of its initial 
value in the first cycle, instead of its low value (7°) obtained without relaxation 
between the Ist and 2nd cycle. 


A NEW METHOD FOR MEASURING AND DEFINING CRIMP 
OF TEXTILE FIBRES* 


M. SHILOH, D. MEJZLER AND E. ALEXANDER 


Institute for Fibres and Forest Products Research, Ministry of Commerce and Industry 


SUMMARY 


The crimp of textile fibres is one of their important characteristics, influencing 
their processing behaviour and the properties of yarns and fabrics made therefrom. 
Despite the importance of crimp, little is actually known about its physical or che- 
mical nature, and methods of measuring and defining crimp are inadequate. 


Existing methods of testing fibres do not take into account the crimp, and there- 
fore may introduce considerable errors in tensile and other measurements. Some 
works deal with the behaviour of the crimped fibres under the application of loads, 
and efforts have been made to express such changes mathematically. No work 
has been done in order to establish a method for measuring the crimp parameters 
in a systematic way so as to be able to distinguish quantitatively between various 
crimps. Such a method is suggested in the present work. 


The statistical nature of the waves has to be taken into account, as the crimped 
fibres are not of a regular shape. It is also necessary to include tests for ascertaining 
the size of representative samples of the fibres, considering the variations in crimp 
within the varieties. 

The main parameter used to describe the geometry of the crimp of a fibre is the 
crimp-diameter. This corresponds to an average amplitude of the fibre in space, 
while it is kept in a static position of stress. 

In order to measure the crimp-diameter, it is necessary to measure the coordinates 
of a sufficient number of points on the fibres. From these measurements a “‘fibre- 
axis” can be defined as the straight line for which the sum of squares of the dis- 
tances of the fibre points from this line is a minimum. The crimp-diameter is then 
defined as: 


D = 2./5/N 


* This work is a part of the research being carried out at this Institute under Grant Number 
FG-Is—103-60 issued by the Agricultural Research Service, United States Department of Agri- 
culture. The experiments performed on wool are a part of a research project carried out under a 
grant of the International Wool Secretariat, London. 

Full details of this work were sent for publication to the “Textile Research Journal”, U.S.A. 


Received May 15, 1961. 
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where S is the mentioned minimum sum and N the number of the points” 
measured along the fibre. The calculation of S according to its definition involves _ 
the solution of third order equations. Therefore the following simplified equation — 


was developed and suggested: 


S = Var(x) + Var(y) — (1/Var(z)) -(Cov7(x,z) + Cov7(y,z)) 


It was found that this equation can, for all practical purposes, serve as a satisfactory | 
approximation for the results obtained when S is calculated according to the de- . 


finition, i.e., from the fibre-axis. 
The crimp-diameter depends upon the size of the intervals, or the number of 


measurements N. When the size of the interval is decreased, N increases, and the | 
measurements become more precise. As a limiting case we may also get the following | 


equation in an integral form: 
lim S/N = lim(S,/N + S,/N) 
and after integrating we get: 


lim(S,/N) = (1/L) Jods — (4/L’) ( pracy +(12/L’) (f zxdz ° fxd _ 


L 
— (12/L’)( fzxdz)’; 
0 

where L is the distance between the fibre endpoints. 

and a similar expression for lim (S,/N). 

These are later introduced for the calculation of the crimp-diameter: 

2 5 ° 
D* = A(lim (S,/N) +im (S,/N)) 
In order to measure the crimp-diameter, a special apparatus had to be built. 


At first a simple uniplanar device was used as a crimp-meter, where the crimp- 
diameters were measured from one plane of the fibres. However, it was found that 


| 
| 


: 
| 


the crimp-diameter varies with the plane of measurement; therefore this method | 
was followed by a rotatory crimp-apparatus: Here the fibre was rotated throug hits _ 


endpoints, and its image was projected on a screen. From this image the “width” 


values can be measured at N intervals, and used for the calculation of D. A series | 
of limitations of the rotatory crimp-apparatus led to the modification of the measuring — 
method and to the development of the method of two perpendicular planes of a 


fibre. Here the fibre’s upper end is clamped to the hook of a torsion balance and 
its lower end to a moving screw, so that the load-elongation data are also easily 
obtained. The coordinates of N fibre points are measured on each screen, when 
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_ the line connecting the endpoints of the fibre is used as the z-axis; S' is then calcu- 
lated and the crimp-diameter is obtained. Also, with the aid of any curve-tracer, 
both perpendicular projections of the crimped fibre can be traced and the x = x(z) 
|} and y = y(z) functions can be obtained and used in the continuous equations. 
The same equations can be used to measure the sums of squares of the deviations 
_ of any curves from their regression lines, either in their continuous or in the non- 
| continuous forms. 

_ This method of measurement and definition was applied to various types of fibres, 
such as cotton and fine wools, and some interesting results have been obtained. 
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i THE STABILITY OF THE MOTION OF A FLUID CONFINED BETWEEN 
1 TWO CONCENTRIC CYLINDERS ROTATING IN OPPOSITE DIRECTIONS* 


LIONEL RINTEL 
Technion-Israel Institute of Technology, Haifa 


ABSTRACT ** 


i 

| G.I. Taylor observed that for cylinders rotating in opposite directions the secondary 
| motion is one of toroidal vortices as in the case of cylinders rotating in the same di- 
| rection. The difference between both cases is that in the case under consideration 
| the perturbation vortices extend only over a part of the gap, limited by the inner 
cylinder and a cylindrical layer of fluid inside the gap. Taylor explains this peculiar 
| behaviour of the perturbation vortices by use of Rayleigh’s nonviscid criterion. 
' According to this criterion the basic solution is always stable when the square of 
circulation increases with the radius. In contradistinction, when the square of cir- 
‘culation decreases with the radius, the basic solution is unstable. This allows to 
| divide the gap into two regions, one always stable, and the second becoming un- 
| stable for secondary motion strong enough to overcome the damping of the pertur- 
bations by the viscosity. This explans qualitatively Taylor’s observations but 
- quantitatively, as Taylor himself remarks, the observed perturbations extend by 
about 30% beyond the boundary indicated by Rayleigh’s criterion. 


* The results reported here are based on part of the work done under the guidance of Prof. 
P. Lieber, University of California, for partial fulfillment of the requirements for a D.Sc. degree 


at the Technion. ; 
** The entire paper is being reviewed for publication in the “Memorial Volume for Prof. E, 


Schwerin” . 
Received May, 15, 1961. 
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This discrepancy can be explained by taking into account the viscosity. Because: 
of the viscosity, tangential stresses will be present on the boundary of the unstable: 
region, as calculated by the criterion. These tangential stresses will exert traction 
on the neighbouring strata of fluid and will bring them into motion. Thus, the un- 
stable region will extend beyond the boundary indicated by Rayleigh’s criterion.. 
The boundary of the extended region will be characterised by the vanishing of the: 
tangential stresses. 

For the determination of the unstable region and checking the theory, the prob-: 
lem of stability of a rotating fluid limited by a rigid inner cylinder and outer’ 
cylinder free from tangential stresses was solved. The extension of the unstable region | 
was determined by comparing the obtained values of the critical Taylor number | 
with those obtained by Chandrasekhar. The obtained values of the critical wave | 


Q , ; : | 

numbers for « = — < —0.50 are in fairly good agreement with those obtained 
if 

by Chandrasekhar. This may be considered as a justification for the separate con- | 


sideration of a part of the fluid with regard to stability. For « < —0.50 a deviation 
from Chandrasekhar’s values for the critical wave number of up to 3% is obtained. 
This deviation is to be explained by the influence of the rigid wall close to the free | 
surface boundary of the unstable region (Wall effect). 


i 
} 
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CHARACTERISTICS OF TURBULENT FLOW WITH SUSPENSIONS 
OF SOLID PARTICLES* 


CHAIM ELATA 
Hydraulics Laboratory, Department of Civil Engineering, 
Technion—Israel Institute of Technology, Haifa 


ABSTRACT 


| The effects of suspensions of nearly buoyant particles on the dynamics of turbulent 
flow in open channels are investigated. 


Results from earlier experiments by Vanoni, Einstein and others indicated, that 
the presence of suspended sediment in turbulent shear flow would change the shape 
| of the velocity profile and the flow resistance factor. These earlier experiments were 
| made with silts and sands. The differences in specific weight between sediment and 
| supporting fluid caused the concentration distribution of the particles to be non- 
| uniform over the vertical. 


This work simplified the correlation of the changes in flow characteristics, by 
| employing neutrally buoyant polystrene particles in the suspension. This elimi- 
nated the buoyant weight of the particles as a variable and caused a uniform 


} concentration of the suspension through the flow. 


Measurements of the resistance factor and velocity profiles in two-dimensional 
steady open channel flow were made, with different concentration of particles. The 
t resistance factors for the flows with suspensions differed only slightly from those 
for clear water flow. It could be found from the measured velocity profile that for 
the flows with suspensions the velocity gradient increased in the outer region of 
| flow, and decreased near the flume bottom, compared with clear water flow. These 


| differences increased with increasing particle concentrations. 
_ To establish the influence of the particle suspensions on the turbulent characteris- 
‘tics of the flow, a turbulence capacitance-gauge was designed. Intensities and spectra 
| of the turbulence component in the mean flow direction could be measured with 
' this instrument. The resulting measurements showed an increase of the intensities 
of the turbulence component, with increasing particle concentrations. The energy 


spectra obtained conformed with the form predicted by Tchen. The results indicated 


. — 


_ * This work forms part of a Research project on the dynamics of open channel flow with suspen- 
| sions of neutrally buoyant particles, by C. Elata and A. T. Ippen, 1960, sponsored by the M.I.T. 
Hydrodynamics Laboratory, Dept. of Civil and Sanitary Eng. 
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a change in the structure of turbulence rather than a “damping” of the turbulence, 
as was previously believed to be the consequence of the effect of the suspension on 
the flow. 

A phenomenological explanation of the observed phenomena was attempted. 
A more recent description of turbulent shear flow by Townsend, rather than Prandtl’ 
Mixing Length Theory, was found to serve best as a background for the observe 
changes in flow characteristics. It was reasoned that the suspended particles coul 
only directly affect the flow in a narrow region near the boundary, where the length: 
scales of the flow are of the same order of magnitude as the particle size. It was: 
assumed that the presence of small rotating particles near the boundary woul 
intensify the production of small scale turbulence, thereby increasing the eddy, 
viscosity. The increased turbulence production near the boundary supposedly changes: 
the equilibrium conditions in the outer region of the flow where, according tos 
Townsend, the momentum transfer is affected by large scale “‘eddies’’. The resul ) 
would be an increase in the velocity gradient in the outer region, which is in agreement 
with the experimental results. 

The volume concentration near the boundary and the particle size and shape: 
appear to be relevant correlation parameters for the changes in flow characteristics 
caused by suspended particles. The effect of the buoyant weight is apparently only 
secondary, causing an increase of the concentration near the boundary. 

The effective viscosity of suspensions of different concentrations was measured. 
The suspensions behaved as Newtonian fluids over the range of measurements. 


SOME EXPERIMENTAL INVESTIGATIONS ON NATURAL 
FREQUENCIES OF RING-FORMED BODIES* 


M. EHRENREICH AND S. BRAUN 
Technion, Israel Institute of Technology, Haifa 


ABSTRACT 

Some methods for measuring the natural frequencies of ring-formed bodies are | 
described. The methods are limited to bodies whose natural frequencies lie in the 
audible range. 

In one method, the ring was excited by an harmonic force, and the frequencies | 
at which the vibrations had maximum amplitude were measured. In another method 
the ring was excited by an impact, and the resulting frequencies were measured. | 
Both methods gave similar results which were in agreement with analytical 
calculations. 


These methods are mostly valuable in case of bodies of irregular shape, like the | 


cage of a roller bearing. Frequently no analytical solution can be found for such a_ 
case, and the natural frequencies can then be measured by the methods outlined. 


* This investigation was a part of research on the noise of ball bearings, in the Faculty of Mecha- 
nical Engineering of the Technion. 


THE VIBRATIONS OF A MEMBRANE IN CONTACT WITH AIR* 


MICHAEL BENTWICH 
Mechanical Engineering Faculty, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


| This paper considers the mechanical system composed of (i) a thin membrane 
» stretched over an aperture in an infinitely wide rigid wall, and (ii) the air in the half- 
| space on the one side of the wall. If all sources of dissipation are negligible, the 
| system is capable of vibrating freely, so that its energy will alternately be potential, 


as compressed air as well as a strained membrane, and kinetic, as a moving mass 
of air or membrane. Solution of the wave-equation, governing air-sound distri- 


( bution together with the membrane-equation, yields the natural frequencies of the 


system. These are expressed as a function of the membrane’s dimensions, the velo- 
cities of wave-propagation in the air and in the membrane, as well as the ratio 


§ of the mass contained in a unit volume of air to that associated with a unit area 


of the membrane. Theoretical analysis, indeed, shows that the membrane’s vibration 


) is associated with a standing sound wave, the intensity of which decreases with 
the distance from the membrane. The solutions obtained reduce to those rep- 
resenting the vibrations of a membrane in vacuo when the above-mentioned 


mass ratio is considered to be extremelly small. 


ON CURVED TUBES WITH VARIABLE CROSS-SECTION 


A. KORNECKI 
Department of Aeronautical Engineering, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 
The achievements in the theory of thin elastic curved tubes and toroidal shells are 


reviewed. 
Underlined are the peculiarities characterising the stress analysis of this type 
of shells and the recen: asymptotical methods of solution of the basic differential 


/ equations. 


Some new problems concerning shells under nonsymmetrical loads, and curved 
tubes with variable cross-section (spiral-casings) are formulated, and methods of 


their solution are suggested. 


* This article is being reviewed for possible publication in the “Journal of the Acoustical Society 
of America’. 
Received May 15, 1961. 
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FRACTURE OF SOLIDS UNDER COMBINED STRESSES 


E. BEN-ZVI 
Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


Problems of safety and fracture of structures are dealt with in “strength of materials” 
over the theoretical background of the classical theory of elasticity. Its basic hypo- 
thesis is the mathematical model of a continuous body, although the “body-model” 
of modern physics is discontinuous. | 

It will be shown here that a continuous body hypothesis is unable to explain some: 
stress-strain relations, especially in cases of solids under combined state of stress. 

This point of view may be proved observing the physical meaning of Poisson’s: 
ratio, where an axially stressed member will be strained also in the cross directions, 
in spite of the non-existence of forces in the lateral direction. Should we even agree; 
to the assumption that a body model will show strains without suitable forces, we 
are unable to explain what is the reason of fracture in the lateral direction under 
direct compression. 

The meaning of fracture is annulation of cohesive forces between particles or} 
molecules of the material. The annulation of cohesion requires the existance of 
internal forces perpendicular to the direction of the applied forces. Such forces: 
are not recognized by the theory of elasticity. Experiments show different phenomena 
in problems of combined stresses!, especially in problems of triaxial compression2, 
where a large gap is observed between theory and experiment. Considering the above 
we may conclude that the theory of elasticity has been formedas an extrapolation of 
the tension test. 

The question is, whether we are in a position to develop a practical suitable theory 
of strength and/or theory of elasticity, based on a discontinuous material body. 
The problem is certainly rather difficult. But as long as the intention is specified, 
we may provide the way in some steps. 

- Observing some structure of real matter, some different forms of discontinuous | 
body could be assumed as proposed by different authors45. E.g., a simple form 
of discontinuous body may be a cubic lattice containing material in spheres. Such | 
representation of matter explains why in simple compression, failure occurs in a 
direction parallel to the applied forces. This model shows clearly that fracture occurs 
due to destruction of cohesive forces, that means by tension only. The value of such 
tension depends on the internal structure of the material in question, which could 
not be defined by the existing theory of elasticity. It seems that such ideas are held 
n “statu-nascendi” by many authors; but they are not yetclearly specified. Future 
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development of this hypothesis may bring us to a more natural method for the 
| solution of engineering problems. 
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SHEAR STRENGTH ANALYSIS OF SPECIMENS HAVING A HIGH 
DIAMETER-HEIGHT RATIO 


(Triaxial Test) 


E. SHKLARSKY AND M. LIvNEH* 
The Highway and Soil Laboratory, Technion-Israel Institute of Technology, Haifa 


ABSTRACT 


| Triaxial Shear tests serve to evaluate the two strength parameters of soil and as- 
| phalts, namely the angle of internal friction (¢) and cohesion (c). Evaluation of these 
| parameters is usually carried out with specimens having a diameter—height ratio 
| of about 0.5, since under these conditions no buckling stresses are developed and 
the influence of friction is kept to a minimum. 

In practice, however, it is not always convenient to prepare specimens with a 
small diameter-height ratio, either because the available material is limited or be- 
cause the thickness of the compacted field layer, from which specimens are cut, 
is small. 

Since strength results of the triaxial and unconfined compression tests run with 
specimens of a high diameter—height ratio are affected by friction developed between 
‘the rough compression plates and the specimens, evaluation of the two parameters (@) 

and (c) by means of such tests is only feasible provided the quantitative effect of 
‘this friction is known. 
_ In this paper, two theoretical solutions are presented, dealing with the strength 
of specimens of high diameter-height ratio, assumed to be compressed between: 
| a) rough plates. b) partially rough plates. 

Laboratory test results show that friction is a function of the roughness of the 
‘compression plates and that “perfect” roughness between specimens and plates 
is unfeasible. 

These two solutions are based on the rupture mechanism theory similar to solu- 
tions common in the plastic theory. 


* This paper is part of the D.Sc. Thesis of the second author under the supervision of the first 
mentioned author. 
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